{"cells": [{"cell_type": "markdown", "id": "88e87d36", "metadata": {"tags": ["module-gfdocean"]}, "source": ["## The mixed layer\n", "[Geophysical Fluid Dynamics of the Oceans](module-gfdocean) \n", "```{index} Mixed layer\n", "```\n", "In this notebook, the mixed layer will be investigated using NetCDF data from the NOAA (National Oceanic and Atmosphereric Administration). "]}, {"cell_type": "code", "execution_count": null, "id": "f6fd1f12", "metadata": {}, "outputs": [], "source": ["import fsspec\n", "import xarray as xr"]}, {"cell_type": "markdown", "id": "2a29da94", "metadata": {}, "source": ["The mixed layer is mixed by the winds on the surface creating a homogeneous layer. It is being mixed by wind stress, eddies, kinematic viscosity and convection. The stronger the winds, the deeper the mixed layer. The mixed layer is deeper in the Winter than in Summer. In the Winter, the surface water cooles so it becomes denser causing convection as it is more dense than the layer below. More turbulence also deepens the mixed layer.\n", "\n", "

\"image

\n", "\n", "\n", "\n", "The mixed layer is between $10$ and $200\\,m$ approximately thick depending on stratification (density gradients), season, wind stress, convection and velocity gradient (turbulence and eddies)."]}, {"cell_type": "markdown", "id": "844089be", "metadata": {}, "source": ["To find out more about this dataset, we can evaluate a few factors: \n", "You can also try this out for yourself to see how much information is contained in this dataset."]}, {"cell_type": "code", "execution_count": 62, "id": "4119f82f", "metadata": {}, "outputs": [{"data": {"text/html": ["
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset>\n", "Dimensions:                       (n_param: 3, n_levels: 77, n_history: 4, n_calib: 1)\n", "Dimensions without coordinates: n_param, n_levels, n_history, n_calib\n", "Data variables: (12/65)\n", "    data_type                     object b'Argo profile    '\n", "    format_version                object b'3.1 '\n", "    handbook_version              object b'1.2 '\n", "    reference_date_time           object b'19500101000000'\n", "    date_creation                 object b'20161104225754'\n", "    date_update                   object b'20170714162347'\n", "    ...                            ...\n", "    parameter                     (n_calib, n_param) object b'PRES           ...\n", "    scientific_calib_equation     (n_calib, n_param) object b'PRES_ADJUSTED (...\n", "    scientific_calib_coefficient  (n_calib, n_param) object b'Surface pressur...\n", "    scientific_calib_comment      (n_calib, n_param) object b'Pressure adjust...\n", "    scientific_calib_date         (n_calib, n_param) object b'20170714162348'...\n", "    crs                           int32 -2147483647\n", "Attributes: (12/49)\n", "    title:                           Argo float vertical profile\n", "    institution:                     CORIOLIS\n", "    source:                          Argo float\n", "    history:                         2018-06-09T01:59:23Z csun convAGDAC.f90 ...\n", "    references:                      http://www.nodc.noaa.gov/argo/\n", "    user_manual_version:             3.1\n", "    ...                              ...\n", "    time_coverage_end:               2015-09-07T12:29:26Z\n", "    time_coverage_duration:          point\n", "    time_coverage_resolution:        point\n", "    gadr_ConventionVersion:          GADR-3.0\n", "    gadr_program:                    convAGDAC.f90\n", "    gadr_programVersion:             1.0
"], "text/plain": ["\n", "Dimensions: (n_param: 3, n_levels: 77, n_history: 4, n_calib: 1)\n", "Dimensions without coordinates: n_param, n_levels, n_history, n_calib\n", "Data variables: (12/65)\n", " data_type object b'Argo profile '\n", " format_version object b'3.1 '\n", " handbook_version object b'1.2 '\n", " reference_date_time object b'19500101000000'\n", " date_creation object b'20161104225754'\n", " date_update object b'20170714162347'\n", " ... ...\n", " parameter (n_calib, n_param) object b'PRES ...\n", " scientific_calib_equation (n_calib, n_param) object b'PRES_ADJUSTED (...\n", " scientific_calib_coefficient (n_calib, n_param) object b'Surface pressur...\n", " scientific_calib_comment (n_calib, n_param) object b'Pressure adjust...\n", " scientific_calib_date (n_calib, n_param) object b'20170714162348'...\n", " crs int32 -2147483647\n", "Attributes: (12/49)\n", " title: Argo float vertical profile\n", " institution: CORIOLIS\n", " source: Argo float\n", " history: 2018-06-09T01:59:23Z csun convAGDAC.f90 ...\n", " references: http://www.nodc.noaa.gov/argo/\n", " user_manual_version: 3.1\n", " ... ...\n", " time_coverage_end: 2015-09-07T12:29:26Z\n", " time_coverage_duration: point\n", " time_coverage_resolution: point\n", " gadr_ConventionVersion: GADR-3.0\n", " gadr_program: convAGDAC.f90\n", " gadr_programVersion: 1.0"]}, "execution_count": 62, "metadata": {}, "output_type": "execute_result"}], "source": ["# evaluate the metadata of this Dataset\n", "ocean_data"]}, {"cell_type": "code", "execution_count": 63, "id": "45d3fa62", "metadata": {}, "outputs": [{"name": "stdout", "output_type": "stream", "text": ["The location of this dataset is at ( -3.6549999999999727 , -10.097 )\n", "The minimum value of temperature in this part of the ocean is 3.397 deg Celcius\n", "The maximum value of temperature in this part of the ocean is 22.737 deg Celcius\n"]}], "source": ["print('The location of this dataset is at (', ocean_data['longitude'].data,',', ocean_data['latitude'].data, ')')\n", "print('The minimum value of temperature in this part of the ocean is', min(ocean_data['temp'].data), 'deg Celcius')\n", "print('The maximum value of temperature in this part of the ocean is', max(ocean_data['temp'].data), 'deg Celcius')"]}, {"cell_type": "code", "execution_count": 64, "id": "4f78639f", "metadata": {}, "outputs": [{"data": {"text/plain": ["Text(0, 0.5, 'Depth from surface in number of levels')"]}, "execution_count": 64, "metadata": {}, "output_type": "execute_result"}, {"data": {"image/png": "iVBORw0KGgoAAAANSUhEUgAAA3gAAAGiCAYAAAC1VGQuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAB4RklEQVR4nO3dd3gc1dXH8e9Rd5Fsy5a7jdyNjXE32DTTew29ExJCElJIJQkBvwQSkkAKNfSe0ItDbzYlgLuNbdx7L5Jlq1j9vn/MyKyFyqqsZnf1+zyPHu3Ozs6cnZ3ds2funTvmnENERERERERiX0LQAYiIiIiIiEjzUIEnIiIiIiISJ1TgiYiIiIiIxAkVeCIiIiIiInFCBZ6IiIiIiEicUIEnIiIiIiISJ1TgiTSBmbkw/tYGHWdLM7NsM5tiZv2DjiUIZnal/95nt9D6Hg/dz/zt78zsykYs6xux++/lMc0SbMNiaZb9yH89U5oprJqWP9mPM6Ha9Kr34TuRWncNsdS0L9S4Dc1srZk93cj1TDGzuL/OUvXtKSISC1TgiTTNxGp/W4F3qk07O7DogpMN3Ay0ygIvCmzB2/feaMRz3/CfuyVk2s1Aixd4xM5+NBkvzmjIqX9g/++cbCKzDR/G20/iXfXtKSIS9ZKCDkAkljnnvgi9b2YlwM7q0+OBmaU650paewyxwN9GjdoHnXM7gB3NG5G0FOfcqhZaz0ZgY0usC8DMEgFzzpW31Dqh5bZnc4ql78lYilUklkTD0UaRuGZm/czsGTPbYWYlZjbfzM6uNs8UvyvXUDN7x8wKzWy9mV3lP36ZmS01swIzm2ZmA6o9f62ZPW1m3zWzlWZWbGZzzezoGuI5ysw+MLN8fz3vmNlB1eaZbmafmtnpZjbPL1x/4D92nZl9bma5ZpZnZl+Y2akhz50MTPPvvhfSVXWy//g3usvV1KXQ7xq10cwmmtlnZrYX+Iv/WBczu9/MNvnbdKmZXRPGe5FmZn83s0X+ttxqZv81s6HV5qvqpnio/97tMbPNZnaXmaVVm7e/mb1hZkX+e/xPILW+WPznnui/tt1+PMvM7KaQxwea2VNmtsbM9prZav91d6pnuXVtz9Fm9okf7wozu7aW157t36/qhve7kPdyipn9wt/2WdWeb36c/6knxibtR7UsM9HMbjWzLf7rm25mw2uZd6SZTTWzXf62/Z+ZHVFtnqptNsnMZvmfq7Vm9qOQeabgtZABlFXFWW11iWZ2ix9Xnr/P9a5n+/zCvM9nSsi0l/zlHxcy7btmVm5mGSExr/VvTyaMbWhmF5rZEn99s83s8Lpiq3rd1V+nv+xbzezH/j6bb2Yf1fQemNnZ/jYv8D9fM83sjGrLus3MbjCzNUApMMJ/LJzvsBPM7M2QfWGRmf3cvEIxdL6LzfuOKzDvc7jQzL4X8nht3Z+/V997amZtzfu85vixvuLvS/V2n7avc8II877zi/x13WIhXYHN6x7szOwcM3vIzHYA20Ie/66ZLfD33Z1m9oiZZVZb10/893+v/3mYbSE5yur/nqqxG6t5n7/pzRmriIRHBZ5IBJlZH2AGMBK4HjgDmAu8FPpjJsQLeF3kzgLmAI+a2R+B7wM3AFcBQ4B/1/Dco4CfAb8DLgRKgLfMbEhIPKcCHwAFwKXAxUA68Ikfa6jBwF3A3cCJ/vPA6/L1MHAecAEwG3jdzE72H58L/NC//WO+7qo6t8aNVLcOwLPAf4CTgX/7P2T/B5wKTPH//xe430J+eNciFe/13uo/7/tAGvCFmXWvYf6ngFXAOcD9/uv6TdWD5v34fg8Y7T92JdAPuLG+F2beOVFTgTV42/EM4G9Au5DZeuK1kvwU7z24BTgWeLO+5dciA2/feRo4E5iFt92+cSAgRFU3vMf5+r18GHgUqMTbJ0OdgLcNHqgnlmyafz+aAvwWeAbvM/Qu3jbej5mNAT4DMoHvAt8CcoD3zWxstdkzgOeAJ/xlTgfuCvmB/jDwiH/78JA4Q/0GGAh8G/iJ//gzdbwOgA+BtsChfsyG9xnfy/7dZY8B5jjn9tSwjHC24RHAz4Hf470PiXjvQ8d64qvNpXifrZ/g7Rt9gdfMbF+PIf9z+jKwHbgCbx94BW+fCHWlv6xf+P83N+A7rL8/37f95z6Bt3/cFhLH4XifhY/w3tvzgIeAcF57OO/pg/7jd+B9hyyrYZ76vAq878f3b7z36aYa5rsbMOAyvO2Gmd0O3Oc//wzgl8BJeHkh0Z/nEuBOvO/YU4BLgBfxPhvhfk81VKNiFZEGcM7pT3/6a6Y/YC3wdMj9R/C6u3WuNt97wPyQ+1MAB1weMq0TUI73wzMjZPqP/XkPqLbeUqBvyLR0IBd4KmTaSuCDarFkADuBf4RMm473431UPa83Aa+r97vAayHTJ/sxHlfDcxwwpdq0bH/6lSHTHvennVlt3t8DxcCgatMf8l9HUgPer0S8H9H5wPUh06/01/1/1eZ/HVgecv+7/nyHVtsmi/3p2XWs+1x/nowGxJuEV0Q4YHS1bbU2zO15dMi0VH+bPVjDa88OmeaAW2uI53F/n7KQaS8DSxv4uWnwflTDMjrh/ej/V7Xpv66+z+H98F8CpFTbF5YAr9awzS6stsz3gHVVr5uvP79J1eareh8+qjb9F/70nvVsk1zgZv/+KLzP5N+Bz0Pm2wLcXse+UOs2xPve2AV0Cpk2zp//4nq29xTAVZvmgBVAcg37+ST/fgbe5+3lepbvgM1Am2rTw/oOq/a4+fvX7/zXmxDyPuTWE0f17RnWe4p3IK4S+FW1+e6i2mezru0L3FBt+kP+9utY7f19pYZ9rwK4qdr0w/z5z/Lv3wPMrSOOer+nqm+jkOnTgek17IuNilV/+tNf+H9qwROJrJPwWlt2m1lS1R/eQCwj/daoUG9V3XDO7cI7wv2F2//o/FL/f/UWty+cc+tDnp/P1wNmYGaDgAHAM9ViKQI+B46stry1zrn51V+QmY01s9fNbBteAVoGHI/3g6a5leMVVaFOwmsVXVPDNu0MDKtrgWZ2vpnNMLM8f/mFQHtqjr/6ICUL8VokqkwENriQcy6dc5XA8/W9MGA+3rZ71szONbOuNcSaYma/Na8L6l5//k/8hxuzvYucc9NCYi3B+0Het/an1Ok+vH3qWD/eHsDp1N96F4n9aAReq0L1bf9stfW2wWsJewGoDNl/DK/1oPrnoAJ4qYZl9gV6hRlbTfsR1LHd/f3oY75urTsG+BLv9Y0zs3QzGwZ0x2vta6zP/e+asGOrx3vOubI6ljcJ7/P2YBjLets5t7fqTkO+w8ysh5k9YGbr8A5+leG13HcEqj5rs4BO5nVvP62BrZb1vaeH4O1TL1Sb78UGrANq3p/bAwdVm/5KtfvH4x0kqL6tZgB7+HpbzQJGmdndZnacmbWttpz51PM91QiNjVVEwqQCTySyugKX4yXI0L+/+o93rjb/rmr3S2uZBl7XwlDb+KZtfP0jtCoxP1JDPKfVEMuWaverupx+gNd950d4P9bGA2/XEE9z2O6cq6g2rStewq/+Gqp+SFV/HfuY2el43e2W4HXtOgQv/h3UHH9utfsl7H9+XQ9q3+51cs6txOt2mYDXFXSrX3geFTLbn/CO5D+N181sAl5XL2qJtz7V9yXwXlOj3jvn3Ey8rpVV5/F9B69Ye6Ku50VoP+rh/6++7avfz8Rrrfs939yHrsP7wR+aG3dVK1hClxlugVfTfgT1v9YPgUP9ovRovPPpZuG1YB/hTyvD67LcWPvF5r4e8KKx70N9r7Xq8xnOAC3Vv4PC+g7z37+p/rRb8Yrj8XzdPTMNwDn3EV63zD54RccOM3vfzA4OI7b6XmfV/ri92nz1fjfUM39t+15t22ol39xWGXz9PjyJ11X9ELyDZLlm9rL55+CG+T3VUI2NVUTCpFE0RSIrB6/F5c+1PL65GdfVrZZpm0JiAe/ckfdrmLe02n1Xwzwn4Z0Xd77zRtEDvMEEGhBnCZBSbVptCbymGHLwfjT9pJbnLKtj3RcCK51zV1ZNMLNk/PNNGmELUNMgHjW9F9/gt6ZNM7NUvO5ItwBvmFm2c26nH++TzrlbQ+Jt38hYI+V+4AEz64VX4L3gnKv+47e65tiPqqv60dgNr4ssIfdD5eF1nbsX78ftN/itZ1U6mVlytSKvapmbiKxpeJ+VI/2/B51z5Wb2CV7R0g+Y6ZwrjHAczWmn/78XsKieeat//sP9DhuA19X0Mufcvuv8+Qd49l+Bcy8CL/qfq8l439Vvm1nvavtBQ1Xtj13xzl+rEtZ3Q7X5V9fw/Or7Xm3b6gRqPrCTA34fW6/F/QHzBm86Ae+cvOfwir5wvqeK+eZ3Onjf6zk1TG9UrCISPhV4IpH1Nl43vsWhXY0i5FAz6+Oc2wBgZul4rT5VXYmW4Z1zM9w5d3sj11H1A3zfj10zG4yX9EOPyFcdzW5TwzLW8c3uRafWMF9t3sZr9VnvnKt+dLw+bfFamEJdhtei0xifA1eZ2aFV3TT91oPzG7IQv9XkQ/9H5mt4P9x3+vFWbz2qPqhJSyil5vcSvMEZ7sAbAKIv8K8wltcc+1F1X+J1tz2f/bssXhg6k3Ou0C+QRuKde1Tfj/hEvEFYQrt6Xgis5+sf2aFx5ocRa7gW4bUu/xKv++nH/vQP8QbD6IPXTbYuDdmGLeEzvHMlr8FrMWqIcL/Datq/kvG2WY2ccwV4g8v0B/6JV5w05XIhM/AKmfPwR//1ndfA5ZwPhL7WC/G2X33F8Xt4BzL6OufeC2dFflfd58zsEOB7NTxe2/fUOqCbmXXxCz7MG+l5CN77XZ8GxyoidVOBJxJZNwEzgY/N7B68Hyed8Aqc/s65bzfjurYB75o3bHsJ3uAS7fAu1ItzzpnZD/FGtEvBO7djJ94R4Ul4BdPf6lnH+3gF0pNmdideN6T/w/uxG9qtbbk/37fNLNePZ5l/XuCzwI1m9ju8a7UdAVzUgNf5d7zR3D4xs7/j/ehrBwwFjnDOnVnHc98GzvKf9zowFm/QmrwGrD/UE3ijm75sZr/Fa1m8Fq9bUZ3MuzzBkXjnaG4AuuC1TGzm6x9vbwNXmNlCvO5L5+C9Vy3tK+BUM3sb7wj7ZufcZgDn3F4zexxvlNiFzrlwftA1x360H+dcnv++/s7M8vEGbBkPXF3D+n+GVyy9Y2aP4LW2dAHGAInOuRtC5s0H/mJmXfDOV7wIOA5vkIyqloiv/P8/N7O3gArn3OwwtkOd/M/sdLyiYJZzbrf/0DS+7uY9rabnhgh7G7YE51y+mf0GuNvMXsIbVTIfbxCZYufc3XU8N9zvsCV4RcdtZlaBV+hdX315ZnaL/9xpeJ+73njfB/Oddz3IprzOZWb2b+AP/kGfOXitrlWtiOG2Dn7Xf/4svK6S38EbMCivnvWvMrM/A/eYN5LyR3gtbX3wznl72Dk3zcwexNv+n+N9fw3GO+j1LoT9PfUCXp55xsz+FjJPVWttncKNNZxliYhH5+CJRJA/6Mk4YAHwR7wjlffjDfLQlIERavIRXteaP+J1r0kDTnbOLQ+J5028ZN0Ob3j3d/COLnfHS/B1cs4txjsKfgDeOS6/witwPq42Xw7e+Uwj/bhm4RVT4J1Xdo//+KvAgXg/KMLi/8idhPeD49f+a3gUb9j/+n4EPIR3Hs4FeJdWOBXvB9fuup5URyyleD9A5uO1pDyB1x3r1jqeVmUB3vvwJ7wfU/f4zz0mpLX3R3jb+Ta89zSdhhXDzeU6vNax/+K9l9WvOVh1/mO9g6tAs+1HNZmCt/9f5i/3BL7+QR263Ll4xV8O3qiG7+K12oyoHgPeIA8X4g3n/xreeW8/cc6Fnmf4Ot77/wO8z9GsOmJsqKp9OvT7Yh5eoV1CPZ/bRmzDiHPO3YNXtPbGK/BewhutcU1dz/OfW+93mP+5PAvYitcN916897V6q98MvBEc/4733fxnvG3UkB4FdbkG77vpV3jn+A3n68tWhPudcybed8xUvMtC3Ip/0K4+zrnf+jEciVcMv4b3nbkL72AFeOdvjsXbf9/DG2n0abz9HcL4nvLP0zsXr9vtq/7r/RnewYWwhBmriISpaohnEYlh5l1k9lPn3KVBxyKtj5ndhndOZE9X8/XYYpLfMnmcc67Oi5KLhMvMfolXSGaHjnpcw3xTgJvxLjlRvVu5iEid1EVTREQaxcxG451n8xO8AUDiprgTaSozOw2vO/58vC6ZR+BdL+/5uoo7EZGmUoEnIiKN9QreOUzv4LU2iMjX8vG6it6A181xE16XYH1WRCSi1EVTREREREQkTmiQFRERERERkTihAk9ERERERCROqMATERERERGJEyrwRERERERE4oQKPBERERERkTihAk9ERERERCROqMATERERERGJEyrwRERERERE4oQKPBERERERkTihAk9ERERERCROqMATERERERGJEyrwRERERERE4oQKPBERERERkTihAk9ERERERCROqMATERERERGJEyrwRERERERE4oQKPBERERERkTihAk9ERERERCROqMATERERERGJEyrwRERERERE4oQKPBERERERkTihAk9ERERERCROqMATERERERGJEyrwRERERERE4oQKPBERERERkTihAk9ERERERCROqMATERERERGJE0lBBxCOLl26uOzs7KDDEBGRCJszZ85O51xW0HHECuVHEZHWI9wcGRMFXnZ2NrNnzw46DBERiTAzWxd0DLFE+VFEpPUIN0eqi6aIiIiIiEicUIEnIiIiIiISJ1TgiYiIiIiIxImYOAevJmVlZWzcuJHi4uKgQ6lVWloavXv3Jjk5OehQRESklYiF/Bg05WcRiWcxW+Bt3LiR9PR0srOzMbOgw/kG5xw5OTls3LiRfv36BR2OiIi0EtGeH4Om/Cwi8S5mu2gWFxfTuXPnqE1eZkbnzp11BFVERFpUtOfHoCk/i0i8i9kCD4j65BXt8YmISHxS/qmbto+IxLOYLvCClJOTw6hRoxg1ahTdu3enV69e++6XlpYGHZ6IiEigbrvtNoYPH87BBx/MqFGjmDFjRq3zXnnllbz44osAfOc73+Grr76qc9n/+te/ePLJJwF4/PHH2bx5c/MFLiIS42L2HLygde7cmfnz5wMwZcoU2rdvzy9+8YtggxIREYkCn3/+Oa+//jpz584lNTWVnTt3hn3w8+GHH653nmuvvXbf7ccff5yDDjqInj17NjpeEZF4EkgLnpmdZGbLzGylmd0QRAwiIiLRKB5y5JYtW+jSpQupqakAdOnShZ49e3LLLbcwfvx4DjroIK655hqcc9947uTJk5k9ezYA7du353e/+x0jR47k0EMPZdu2bYB3YPWOO+7gxRdfZPbs2VxyySWMGjWKN954g7PPPnvfst577z3OOeecFnjFIiLRo8Vb8MwsEbgXOB7YCMwys6nOubr7Y9Th//67mK8272muEAEY1jODm08f3qzLFBERqUtz58ig8uMJJ5zALbfcwuDBgznuuOO44IILOOqoo7juuuu46aabALjssst4/fXXOf3002tdTmFhIYceeii33XYbv/rVr3jooYe48cYb9z1+7rnncs8993DHHXcwbtw4nHP8/Oc/Z8eOHWRlZfHYY49x1VVXNc8LFxGJEUG04E0AVjrnVjvnSoFngTMDiENERJrJ4s27+WJ1TtBhxIO4yJHt27dnzpw5PPjgg2RlZXHBBRfw+OOPM23aNA455BBGjBjBhx9+yOLFi+tcTkpKCqeddhoAY8eOZe3atXXOb2ZcdtllPP300+Tl5fH5559z8sknN9fLEpFWYOaaXBZt2h10GE0SxDl4vYANIfc3AodUn8nMrgGuAejbt2+dC1RLm4hIsB76eDXzNuTx0S+PDjqUWFdvjoyV/JiYmMjkyZOZPHkyI0aM4IEHHuDLL79k9uzZ9OnThylTptR7qYLk5OR9I14mJiZSXl5e73qvuuoqTj/9dNLS0jjvvPNIStJwAyISvpteW0TfzLY8ePm4oENptCBa8Goam/gbnfCdcw8658Y558ZlZWW1QFgiItJYZRWO5EQNzNwM6s2RsZAfly1bxooVK/bdnz9/PkOGDAG88/EKCgr2jZrZVOnp6eTn5++737NnT3r27Mmtt97KlVde2SzrEJHWo6i0grYpiUGH0SRBHNbaCPQJud8b0PjGIiIxrLSikqQEXVusGcRFjiwoKOBHP/oReXl5JCUlMXDgQB588EE6duzIiBEjyM7OZvz48c2yriuvvJJrr72WNm3a8Pnnn9OmTRsuueQSduzYwbBhw5plHSLSehSVVtA2NbZb/q2mEawiukKzJGA5cCywCZgFXOycq7Uj/rhx41zViFpVlixZwoEHHhjJUJtFrMQpItIUVz42k9zCUqZed3iTlmNmc5xzsdsvpokamiNjOT9G0nXXXcfo0aO5+uqra51H20lEajLspre5eEJfbjwt+g4QhZsjW7w8dc6Vm9l1wDtAIvBoXcWdiIhEv/IKpxa8ZqAc2XRjx46lXbt23HnnnUGHIiIxqLS8kpSk2D7lIJD2R+fcm8CbQaxbRESaX3llJUkJsZ0Qo4VyZNPMmTMn6BBEJIa1SU5kb1lF0GE0ibKxiIg0WWUlJKoFT0REYlzb1ET2lqrAC0xLnz/YUNEen4hIcymvrFSBF0WUf+qm7SMitWmbkkShCrxgpKWlkZOTE7Vf0s45cnJySEtLCzoUEZGIq6h0KvCiRLTnx6ApP4tIXVKTEiiO8S6aMTsGaO/evdm4cSM7duwIOpRapaWl0bt376DDEBGJuNIKF/MnpceLWMiPQVN+FpHamMX+wcqYLfCSk5Pp169f0GGIiAhQVlFJii50HhWUH0VEWjdlYxERabKyikqSEmP/qKeIiEisU4EnIiJNVl6hc/BERCT2lVdUxvx1XVXgiYhIk1U6R2IcnLcgIiKtW1FpBW1TYvYsNkAFnoiININKpxY8ERGJfUWl5bRNSQw6jCZRgSciIk1W6eJj5DEREWndikoraJuqAk9ERFo55xyq70REJNY5B0ZsJzQVeCIi0mSVDtRDU0REYl2blET2lpYHHUaTqMATEZEmc87F/BFPERGRdimJFJZWBB1Gk6jAExGRJnOgLpoiIhLz2qYmUaQWPBERae0qKx0JqvBERCTGJScmUFbhgg6jSVTgiYhIkzmHCjwREZEooAJPRESarNI5DbIiIiISBVTgiYhIk1U4R4IqPBERkcCpwBMRkWah8k5ERCR4KvBERERERETihAo8ERERERGROKECT0REREREBHAuti+RACrwREREREREACgtryQtOTHoMJpEBZ6IiIiIiAhQWFpOuxQVeCIiIiIiIjGvqLSCNirwRESkNcsrKqWsQtfBExGR2FdR6Vi2NZ/isoqgQ2k0FXgiItIkf3lnGQCnH9wz4EhERESa5pcnDuGzVTlc9sgM8opKgw6nUVTgiYhIo81bv4v/zFzPlZOyGdYzI+hwREREmuSqw/px90WjWbBhN9+6/zM25BYFHVKDqcATEZFGqah03PjqIrqmp/LT4wYFHY6IiEizOH1kT566egI78ks4+77PmLt+V0xdPiEp6ABERCT2FJdVcMNLX7J48x7uvmg06WnJQYckIiLSbA7p35mXvj+JKx+bxTn3fUaPDmmMz85kfL9Mxmd3YnDX9Kg99zxiBZ6ZPQqcBmx3zh3kT8sEngOygbXA+c65XZGKQUREmt+O/BKueWo289bn8csTh3DawT2CDinmKEeKiES/Qd3SmXrdYfx3wWZmrdvFF6tzmLpgMwAd2iQz7oBOjMvOZEK/Tozo1ZGUpOjoHGmRam40syOBAuDJkOT1FyDXOXe7md0AdHLO/bq+ZY0bN87Nnj07InGKiEj4vtq8h+88MYtdRWX8/YKRnHRQ8xZ3ZjbHOTeuWRcahZorRyo/ioi0HOccG3L3MnNtLrPX5jJzbS6rdxQCkJqUwMg+HZngt/KN6dux2Xu3hJsjI9aC55z72Myyq00+E5js334CmA7UW+CJiEgwyisqyS0qJbewlIUbd3Pz1MVkpCXzwrUTOahXh6DDi1nKkSIiscfM6Nu5LX07t+Xcsb0B2FlQwuy1u5i1NpdZa3O5/6NV3DNtJQkGB/bIYHx2Jof2z+SEYd1brEtnS5+D1805twXAObfFzLrWNqOZXQNcA9C3b98WCk9EJL6VlleSW1hKTmEJuYVe4ZZT8PX9nILSr6cXlrJ7b9l+zx/ZpyMPXTaWrhlpAb2CuBZWjlR+FBEJRml5JTmFJezML2VnQYn/9/XtjLRkenVsw/rcIiodLN68h8Wb9/D4Z2t59MpxHDO0W4vEGbWDrDjnHgQeBK8LSsDhiIhEpeKyCnIKS8n1i7SqAi2nsJTcqqLNL9hyC0rJLymvcTkJBpntUvb9HdgjY9/tzu1T6Nwulc7tUxjdtyOpSYkt/CollPKjiEjzKSotZ2d+KTsKSsipVrDlFHjTdxaUsDO/hD3FNefQNsmJdElPoUv7VAZ3S+ewgZ3p0j6VLu293NmjQxpj+nZqsdfU0gXeNjPr4R+Z7AFsb+H1i4hELeccRaUV+7eoFZbu19JWvWgrKq2ocVnJiUantl8XaH06dfRut0shs73/v13qvmkd2iRH7WhgrYhypIhIEznn2LO3fF9hllMQ2tq2fwG3M7+UvWU159GMtCS6pKfSpV0qQ7un02Vgl30FW1XxluXfb5caXW1mLR3NVOAK4Hb//2stvH4RkRZXWl7J5ry9rM8tYsvuveys1g0yt7DEb4ErpaS8ssZlpCYl7CvOMtul0j+r/dctbCEtbVVFW0ZaEmYq2GKMcqSISC0KSsrZtqeY7XtK2J5fzI78Er/Vbf+CLaewhLKKb3ZuMIPO7bweKV3SUxjTt9N+BVtWSItb5/YpMd1bJZKXSfgP3sniXcxsI3AzXtJ63syuBtYD50Vq/SIiLcU5x478EjbsKmJ9bhEbcveyIde7vXHXXrbs3ktltVzTNiVxX3GW1T6VId0y/AKtWtHWLpXM9im0S0lUwRZHlCNFRDx7SyvYtqfY+8svYfueYrbnl+ybtn2Pd7uwhh4ryYm2r2Dr0j6Vod29XFpVrIUWcJntUkhsJT1VIjmK5kW1PHRspNYpIhIpBSXl+4q2DVV/u/b6RVwRxWX7t7x1y0ilT6e2TOiXSZ/MtvTp1Ia+mW3p1akNXdqnkpYcu0cGpemUI0Uk3hWXVbBjX6Hm/88vZseeErblfz0tv4bz2lKSEuiWkUq39DQO7JHBUUOy6JaRtm9a14xUstqnkdFGvVVqEl0dRkVEAlJWUcmWvGKvgNv1zUIut7B0v/nTU5PondmWAVntmDw4i76d29KnU1v6ZLald6c2KuBERCQulZZXsqOgxG9dCyne/K6T2/0CLq+o7BvPTU40uqZ7hdqgru05bEBnumakfV28ZaTRLV2FW1OpwBORVqWsopL5G/KYuSaXdTmFbMj9+ty40G6USQlG705t6JPZlpN6dfCLN68Vrk+ntnRsm6zkIyIicWnL7r3MWruLFdvy92uB255f8o0DngCJCUbX9FS6ZqRxQGev90q3DO9+13S/cMtIo5NyZ4tQgScicW9DbhEfLd/Bx8t38PmqnH2XCshKT6VPpzaMz+5En8xefldK7wKm3TPSWk1ffRERab0qKx0rthcwa20us9fmMmvtLjbl7QW8S+hk+QVa705tGHNAJ7qlf93a1tX/n9k2RSMxRxEVeCISdwpKyvl8VQ6frPCKurU5RQD06tiG00b24MhBWUwa0IUObZMDjlRERKRlFZdVsHDTbr+g28Xstbn7ru+WlZ7K+OxOXH14P8ZnZzK0RzrJiQkBRywNpQJPRGJeZaVj8eY9fLxiBx8t38Hcdbsor3S0SU5k4oDOXDkpmyMGZ9G/Szt1DRERkVYlr6iUOet2Mcsv5r7cuJvSCm9gsAFZ7ThlRA/GZWcyPrsTfTPbKk/GARV4IhJTissqWJdTxJqdhazNKWTx5j38b+XOfecEDO+ZwXeP7M8Rg7ow9oBOMX0dGxERkXAUl1WwdXcxW3Z7lxbYstsbNGzuul0s25YPeAOcHNSrA1cels24AzoxLjuTzHYpAUcukaACT0SiTnFZBetz/SLOL+TW7ixibU4hW3YX7zdvt4xUJg/O4sjBWRw+qAtd2qcGFLWIiEjzcs6xp7h8X9G2dfdetu4uYeuevf79YrbuqXnEyo5tkxnZuyOnj/Ra6Eb27kibFB30bA1U4IlIIIrLKtiQ+3VL3NqcIq+Y21nIlj3FuJARLTPbpZDduS0TB3Qmu3M7sru0o1/ndhzQpS0ZaTqPTkREYk9lpSOnsHRfkbZ1t1+07fELN/92UQ0X+O7SPpXuHVLp3akt47Mz6d4hje4ZafTokEY3/3a7VP3Mb63qfefN7DzgbedcvpndCIwBbnXOzY14dCIS00rKq4q4ItblFH5dzO0sYvPuvfsVcZ3aJpPdpR2H9K8q4trSr0s7Dujcjg5tVMRJdFKOFJGalFVUsj2/ZF+L25bde0MKOa81bnt+MWUVbr/nJSUY3TLS6N7Bu8D30UO7ekWbX7x175BG1/Q0UpI08InULpzS/vfOuRfM7HDgROAO4H7gkIhGJiIxJa+olFlrdzFrbS5fbd7D2pxCNuftf225Dm28Im58dieyu/TeV8D169xOI1pKrFKOFGmlSssrWbWjgGVb81m2LZ9V2wvY6nel3FlQst9BTIC05AR6dGhD94w0DumXSbcOftHmF3TdO6TRpV2qLjcgTRZOgVfVLnwqcL9z7jUzmxK5kEQkFmzZvZeZa3KZtTaXWWu+Pok7JTGBoT3SGXtAJ84Z05t+XdqS3bkd/bq0o2NbncwtcUc5UiTOVVY6NuXtZenWfJZt3cPSrfks35bP6h2FlPtHMZMTjezO7ejZsQ3DemTs1+LWvUMaPTLakNEmSSNUSosIp8DbZGYPAMcBfzazVEDtwiKtiHOOVTsK/WIul5lrc9m4y7sIaruURMZmZ3L6yB6Mz85kZJ+OpCXrJG5pNZQjReJITkHJvha5ZVvzWbo1nxXb8ikMOQ+ud6c2DO2ezvHDujGkewZDuqXTr0s7dZuUqBFOgXc+cBJwh3Muz8x6AL+MbFgiEqTyikq+2rJnXwvd7LW7yPEvQ9C5XQrjszO56rB+TMjO5MAe6STpIqjSeilHisSgotJyVmwr2FfELd/m/d9ZULJvnk5tkxnSPZ3zxvVhSPd0hnRPZ3C3dNpr8BKJcrXuoWaWGXJ3esi0EmB2ZMMSkZZWXFbBu19t45W5G5m5Jnff0co+mW04akgWE7IzGd8vUxcLF0E5UiSWOOdYl1PE56tz+GxVDl9uzGN9btG+c+TSkhMY3C2dyUOyGOoXckO6p5PVPlX5TmJSXYcg5gAOqGnPdkD/iEQkIi3GOceXG3fzwpwNTJ2/mT3F5fTskMY5Y3ozvl8mE/yhl0XkG5QjRaLYpry9fL4qx//byWb/GqpZ6amMO6ATZ4/u5RdzGfTNbEuiBjaROFJrgeec69eSgYhIy9mRX8Kr8zbxwpwNLN9WQGpSAicd1J3zxvZh0oDOGsFLpB7KkSLRZXt+MZ+vyuELv5VuXU4R4HWznDigM9/v35mJA7owIEu9UCT+hXMdPAMuAfo55/5gZn2B7s65mRGPTkSaTVlFJR8u3c4Lszcyfdl2yisdo/p05LazD+K0g3vqWnMijaAcKRKMXYWlzFjjFXOfr8phxfYCANJTkzikf2cun5jNpAGdGdItXQctpdUJ5yzR+4BK4BjgD0A+8BIwPoJxiUgzcM6xYONuXp23if8u2ExOYSlZ6alcfXg/zh3bm0Hd0oMOUSTWKUeKtICyikrmrtvFR8t38NHyHXy1ZQ/OQZvkRMb3y+RbY3szaUBnhvfsoO6W0uqFU+Ad4pwbY2bzAJxzu8xMF7MSiWLrcgp5dd5mXp2/iTU7C0lJSuDYoV05d2xvjhqcpVEvRZqPcqRIhGzO2+sVdMt28L+VO8kvKScxwRh7QCeuP24wkwZ05uDeHXV5ApFqwinwyswsEe+kccwsC+9opYhEkdzCUt74cjOvzNvE3PV5mMEh/TK59qj+nHRQD3XBFIkM5UiRZlJSXsHstbuYvmw7Hy3fwfJtXrfLHh3SOG1kD44anMWkgV3ISFM+E6lLOAXeXcArQFczuw04F7gxolGJSFj2llbw/pJtvDpvEx8t30F5pWNIt3RuOHkoZ4zsSc+ObYIOUSTeKUeKNMH6nCI+Wr6d6ct28NmqHPaWVZCSmMCEfpmcN7YPRw3JYlDX9hoYRaQB6i3wnHPPmNkc4Fi84aDPcs4tiXhkIlKrnQUlPPrpGp76fB35JeV0y/DOqztrdC8O7JERdHgirYZypEjDOOdYujWftxZu4c1FW1npD47SN7Mt543zTiM4tH9n2uli4iKNFs4omv8EnnPO3dsC8YhIHTbn7eXBj1fz7Kz1lJRXcvJB3bn0kAM4pH9nnVQuEgDlSJH6OedYvHkPby7cwluLtrJmZyEJBof068wlh/Rl8pCu9OvSLugwReJGOIdH5gI3mtlgvG4ozznnZkc2LBEJtXZnIfdPX8XL8zbiHJw1uhffnzyAAVntgw5NpLVTjhSpgXOOLzfu5s1FW3hr4VbW5xaRmGBMGtCZ7x7RnxOGd6NL+9SgwxSJS+F00XwCeMLMMoFvAX82s77OuUERj06klVu6dQ/3TVvF619uJikxgYsm9OWaI/vTu1PboEMTEZQjRUJVVjrmbcjjLb+lblPeXpISjMMGduG6owdy/LBudGqnQWZFIq0hHZwHAkOBbOCriEQjIoB3jt2f31rKC3M20i4lke8e0Z+rj+hH1/S0oEMTkZopR0qrlVdUyotzNvL0F+tYm1NESmICRwzqwvXHD+b4A7vRoa1GvRRpSeGcg/dn4BxgFfAc8AfnXF6E4xJplSoqHc/MWMcd7yyjqLSC7x3Zn+9PHkDHtjriKRKNlCOlNVu4cTdPfbGW1+ZvpqS8knEHdOJHxwzi+OHddCkDkQCF04K3BpjonNsZ6WBEWrM563Zx02uLWLx5D4cN7Mz/nTGcgV3Tgw5LROqmHCmtSnFZBW98uYUnv1jHgg15tElO5Jwxvbns0AMY1lOjOItEg3AKvAeBi82sv3PuFjPrC3R3zs2McGwircLOghL+8vZSnp+9ke4Zadxz8WhOHdFD1/wRiQ3KkdIqrMsp5N8z1/P8rA3sKipjQFY7ppw+jHPG9lZrnUiUCafAuxeoBI4BbgHygZeA8RGMSySuOef4YnUuz85az1uLtlJZ6fjeUf358TGDdO0fkdiiHClxqayiktlrdzFt2XamLd3Oiu0FJCYYxx/YjcsnHsDEAZ11IFIkSoXzS/IQ59wYM5sH4JzbZWb1nhBkZn2AJ4HueMnvQefcP/2Rxp7DOxF9LXC+c25XI+MXiSk78kt4ae5Gnpu1gTU7C0lPS+LC8X24YlK2LnkgEpsanCOVHyVabd9TzPRlO5i2bDufrthJfkk5yYnGhH6ZXDC+D6ce3IMeHdoEHaaI1COcAq/MzBIBB2BmWXgJqT7lwM+dc3PNLB2YY2bvAVcCHzjnbjezG4AbgF83KnqRGFBZ6fhk5U6enbme977aRnmlY0J2Jj86ZiCnjOhBWnJi0CGKSOM1JkcqP0pUqKh0zN+Qx/Rl25m2bDuLNu0BoHtGGqeN7MHkIV05bGAX2qtniUhMCecTexfexVu7mtltwLnAjfU9yTm3Bdji3843syVAL+BMYLI/2xPAdJTAJA7tyC/hPzPX89ysDWzK20tmuxSuOiybC8b3ZWBXtdaJxIkG50jlRwlSYUk57y/ZxrSl2/lo+Q52FZWRYDD2gE788sQhHD2kKwf2SFf3S5EYFs6Fzp8xsznAsYABZznnljRkJWaWDYwGZgDd/OSGc26LmXWt5TnXANcA9O3btyGrEwmUc46pCzZz02uL2b23jMMHduE3pwzl+GHdSE1Sa51IPGlqjlR+lJayI7+EJz5by1NfrGP33jIy26Vw9JCuHD20K0cM6qLL8YjEkVoLPP9cgCrbgf+EPuacyw1nBWbWHu+E85865/aEe0TIOfcg3uhkjBs3zoX1JJGA5RSUcOOri3hr0VbG9O3IX849WJc6EIlDzZEjlR+lJazeUcBDn6zhpbkbKauo5IRh3bj68P6MO6ATCQlqpROJR3W14M3BO6eg6tNflUTMv92/voWbWTJe8nrGOfeyP3mbmfXwj072wEuMIjHv7UVb+d0rC8kvLueGk4fy3SP6k6jkKRKvmpQjlR8l0uau38UDH63i3a+2kZyYwLfG9Oa7R/Sjvwb0Eol7tRZ4zrl+TVmweYciHwGWOOf+FvLQVOAK4Hb//2tNWY9I0PKLy7jptcW8Mm8Tw3tm8O/vjmJId7XaicSzpuRI5UeJpOnLtnPvtJXMWruLDm2S+eHkgVwxKZus9NSgQxORFhLJYZEOAy4DFprZfH/ab/ES1/NmdjWwHjgvgjGIRNSe4jIue2Qmizbt5ifHDuK6YwaSnJgQdFgiEt2UH6XZFZaU83//XczzszfSq2MbbjptGBeM76Nrq4q0QhH71DvnPuXrrivVHRup9Yq0lPziMi5/ZCZfbd7N/ZeM4YTh3YMOSURigPKjNLf5G/L46bPzWJdbxHVHD+Qnxw3SwUaRVqyuQVb6OefWtGQwIrEiv7iMyx/1Wu7uU3En0uooR0o0qKh03D99JX9/fwXdM9J49ruHckj/zkGHJSIBq+vwzosAZvZBC8UiEhMKS8q54tGZLNy4m3suVnEn0kopR0qgtu8p5qIHv+COd5dzyogevPmTI1TciQhQdxfNBDO7GRhsZj+r/mC1E8NFWo2HPlnN3PV53H/JGE46SMWdSCulHCmBcc7x8xcWsHDTbv5+wUjOGtVLFyYXkX3qasG7ECjGKwLTa/gTaXXyi8t49NM1HD+sGyeP6BF0OCISHOVICcyr8zfxyYqd3HDyUM4e3VvFnYjsp67LJCwD/mxmXzrn3mrBmESi1pOfr2NPcTk/PmZQ0KGISICUIyUouYWl/OH1JYzu25FLDz0g6HBEJAqFM8TSZ2b2NzOb7f/daWYdIh6ZSJTZvbeMhz9ZzdFDshjRWx8BEQGUI6UFlVdU8vvXFpFfXMbt5xxMYoJa7kTkm8Ip8B4F8oHz/b89wGORDEok2uzeW8blj8ygoKScnx43OOhwRCR6KEdKi9ict5eLH5rBG19u4cfHDGJId/UEFpGahXMdvAHOuW+F3P+/kAuzisS9vKJSLntkJsu25nP/JWMZ2adj0CGJSPRQjpSIe++rbfzyxQWUllfyt/NHcs6Y3kGHJCJRLJwCb6+ZHe5fmBUzOwzYG9mwRKJDbmEplz48g5U7CnjgsrEcPbRr0CGJSHRRjpSIKSmv4E9vLuXxz9YyvGcGd180mv5Z7YMOS0SiXDgF3rXAkyHnFOwCrohcSCLRIbewlIsf+oI1Owt5+PJxHDk4K+iQRCT6KEdKRKzLKeQHz8xl8eY9fPuwfvz65CGkJiUGHZaIxIB6Czzn3AJgpJll+Pf3RDwqkYDt3lvGZY/MYM3OQh69cjyHDewSdEgiEoWUIyUSpi3dzk+enYeZ8fDl4zhuWLegQxKRGBJOCx6gpCWtR0FJOVc+NpPl2/J56PJxKu5EpF7KkdIcKisdd324gn9+sIIDu2fwwGVj6ZPZNuiwRCTGhF3gibQGe0sr+M4Ts/hy427uu2QMk4fonDsREYm83UVlXP/8fD5cup1zxvTij2ePIC1ZXTJFpOHqLPDMLAE41Dn3WQvFIxKY8opKvv/MHGasyeUfF4zixOHdgw5JRKKYcqQ0lw25RVz6yAw27drLH84czqWHHoCZrnEnIo1T53XwnHOVwJ0tFItIoG59YwnTl+3gtrNGcOaoXkGHIyJRTjlSmsPW3cVc8vAM8orKePaaQ7lsYraKOxFpknAudP6umX3L9G0jceyZGet4/LO1fOfwflx8SN+gwxGR2KEcKY2WU1DCJQ9/QW5hKU9+ewLjsjODDklE4kA45+D9DGgHVJjZXsAA55zLiGhkIi3k81U53PzaYiYPyeI3pxwYdDgiEluUI6VRdheVcdkjM9mUt5cnv30II/t0DDokEYkT4VwmIb0lAhEJQk5BCT94Zg7ZXdpx10WjSUzQQXgRCZ9ypDRGZaXju0/OZuX2Ah6+YhwT+qnlTkSaT71dNM1zqZn93r/fx8wmRD40kch7Yc5GdhWVcc/Fo8lISw46HBGJMcqR0hhTF2xm5tpcbj37II4cnBV0OCISZ8I5B+8+YCJwsX+/ALg3YhGJtJDKSsezM9czoV8mQ7urN5WINIpypDRIcVkFf31nGcN7ZnDumN5BhyMicSicAu8Q59wPgWIA59wuICWiUYm0gC9W57A2p4iLJ2hQFRFpNOVIaZCnv1jHpry9/ObkA0nQaQEiEgHhFHhlZpYIOAAzywIqIxqVSAt4ZsZ6OrZN5qSDdL07EWk05UgJ267CUu7+cCVHDs7i8EFdgg5HROJUOAXeXcArQDczuw34FPhjRKMSibANuUW8tWgLF4zrQ1pyYtDhiEjsUo6UsN353jIKSsr57SlDgw5FROJYOKNoPmNmc4Bj/UlnOeeWRDYskch65NM1JJhx1WH9gg5FRGKYcqSEa/Hm3fx7xnoun5it875FJKLCuQ4eQFugqgtKm8iFIxJ5eUWlPD97A2eM6kn3DmlBhyMisU85UurknGPK1MV0bJvC9ccNDjocEYlz4Vwm4SbgCSAT6AI8ZmY3RjowkUj55wcrKCqt4Joj+wcdiojEOOVICccjn65h1tpd/PqkIXRoq0vyiEhkhdOCdxEw2jlXDGBmtwNzgVsjGZhIJLw0ZyOP/W8tV0w8QF1kRKQ5KEdKnV6bv4lb31jCicO7cd7YPkGHIyKtQDiDrKwFQvuxpQKrIhKNSATNXb+L37y8kEkDOnPjacOCDkdE4sNalCOlFh8v38EvXljAIf0y+eeFo3VZBBFpEbW24JnZ3XjnE5QAi83sPf/+8XijhInEjK27i/neU3Po3iGNey8eQ3JiOMc2RERqphwp9VmwIY9rn57DwK7pPHTFOI3YLCItpq4umrP9/3PwhoCuMj1i0YhEyM1TF1FUUs4z3zmETu10DWIRaTLlSKlVWUUl3396Dp3bp/DEVePJSNN5dyLScmot8JxzTzRlwWaWBnyM110lCXjROXezmWUCzwHZeF1bznfO7WrKukTqsqe4jGlLd3DZxAMY3C096HBEJA40JUcqP8a/D5ZsY/PuYh6+fBxdMzRas4i0rHBG0TzNzOaZWa6Z7TGzfDPbE8ayS4BjnHMjgVHASWZ2KHAD8IFzbhDwgX9fJGI+XLKd0opKThnRI+hQRCTONDJHKj/Guae/WE/PDmkcPbRr0KGISCsUzolI/wCuADo75zKcc+nOuXqHH3SeAv9usv/ngDPxhpTG/39WQ4MWaYg3Fm6he0Yao/t0DDoUEYk//6CBOVL5Mb6t2VnIpyt3ctGEviRqUBURCUA4Bd4GYJFzzjV04WaWaGbzge3Ae865GUA359wWAP+/Dm9JxExftp1pS7dz6sE9NHqZiERCo3Kk8mN8yi8u41cvLiA50bhgvC6JICLBCOc6eL8C3jSzj/C6lQDgnPtbfU90zlUAo8ysI/CKmR0UbmBmdg1wDUDfvn3DfZrIPl+szuF7T81hSPd0fnzsoKDDEZH41KgcqfwYf/KKSrn80Zl8tXkP/7xwtM69E5HAhNOCdxtQhHedn/SQv7A55/LwRhY7CdhmZj0A/P/ba3nOg865cc65cVlZWQ1ZnQjzN+Rx9eOz6JPZlie/PYEObTSCmYhERJNypPJjfNiRX8KFD37B0q35PHDZWE49WOd8i0hwwmnBy3TOndDQBZtZFlDmnMszszbAccCfgal45yvc7v9/raHLFqnL0q17uOLRmWS2T+Hpqw+hc/vUoEMSkfjV4Byp/Bhftu4u5uKHv2BLXjGPXjGewwd1CTokEWnlwinw3jezE5xz7zZw2T2AJ8wsEa+l8Hnn3Otm9jnwvJldDawHzmvgckXqdNNri0lJSuDf3zmU7h3URUZEIqoxOVL5MY5MmbqYrbuLefLqCYzPzgw6HBGRsAq8HwK/MrMSoAwwvEHA6hsl7EtgdA3Tc4BjGxGrSL227N7LzDW5/Oz4wfTJbBt0OCIS/xqcI5Uf48fanYW889VWfjh5oIo7EYka9RZ4zjldGVpixhtfbgHg9JE9A45ERFoD5cjW7ZFP15CckMDlkw4IOhQRkX3qLfDM7MiapjvnPm7+cESa5r8LNjOiVwf6dWkXdCgi0gooR7ZeuYWlvDBnA2eP7kXXdJ0OICLRI5wumr8MuZ0GTADmAMdEJCKRRvpo+Q4WbNzN708bFnQoItJ6KEe2Und9sIKS8kq+e2S/oEMREdlPOF00Tw+9b2Z9gL9ELCKRRthbWsGNry6kf1Y7Lj1U14USkZahHNk6Lduaz1NfrOPiCX0Z2FW9dEUkuoTTglfdRiDsC7KKtIS7P1zBhty9/Oe7h5KalBh0OCLSeilHxjnnHDdPXUR6WhK/OGFI0OGIiHxDOOfg3Q04/24CMApYEMGYRBpkxbZ8Hvx4NeeO7c3EAZ2DDkdEWhHlyNbnzYVb+WJ1LreedRCd2qUEHY6IyDeE04I3O+R2OfAf59z/IhSPSIP9Z+YGEhOM355yYNChiEjroxzZyjz4yWoGdW3PRRN0OoCIRKdwzsF7oiUCEWmsz1btZFx2JzJ1JFVEWphyZOuyakcBCzbkceOpB5KYYEGHIyJSo4T6ZjCzw8zsPTNbbmarzWyNma1uieBE6rOzoISlW/OZNKBL0KGISCukHNm6vDJ3EwkGZ+haqyISxcLpovkIcD3esM8VkQ1HpGE+W5UDwGEDVeCJSCCUI1uJikrHK/M2cdjALnTN0HXvRCR6hVPg7XbOvRXxSEQaqLyikvumraRXxzYc1DMj6HBEpHVSjmwlnp21nk15e/n9aTrfW0SiWzgF3jQz+yvwMlBSNdE5NzdiUYmE4d8z17N0az73XTKGpMR6exuLiESCcmQrkFdUyh3vLOPQ/pmcOLx70OGIiNQpnALvEP//uJBpDjim+cMRCU9OQQl3vLOMSQM6c/JBSrYiEhjlyFbgzneXs6e4nClnDMdMg6uISHQLZxTNo1siEJGG+Pv7yykqreD/lGxFJEDKkfFv2dZ8npmxjssnZjO0u04HEJHop35tEpPe+2obJx3UnUHd0oMORURE4tiLc7xrrf70uEFBhyIiEhYVeBJztu0pZtueEsb07RR0KCIiEsecc7yzeBuTBnShY1tda1VEYoMKPIk5CzfuBmBE7w4BRyIiIvFs6dZ81ucWcZLO9RaRGBLOICuY2SQgO3R+59yTEYpJpE6frtxJgsGwHjoXQkSCpxwZv56btQEzOO7AbkGHIiIStnoLPDN7ChgAzOfri7g6QMlLWtxzs9bz+GdrOWdML9qlhnV8QkQkYpQj49fDn6zm8c/WcuH4PmSlpwYdjohI2ML5hTwOGOacc5EORqQuby/awm9eXshRg7O4/ZyDgw5HRASUI+PSf2au59Y3lnDyQd259ayDgg5HRKRBwjkHbxGgzucSqE9X7OTH/5nP6L6duP/SMaQk6fRREYkKypFx5rX5m/jtKwuZPCSLf144mqRE5RsRiS3htOB1Ab4ys5lASdVE59wZEYtKJMSmvL1c89Rs+me149ErxtM2RV0zRSRqKEfGkQUb8vjZ8ws4pF8m/7p0rA4mikhMCueX8pRIByFSlyc+W0tJeSUPXT6ODm2Tgw5HRCTUlKADkObz7xnrSUtK4KHLx5GWnBh0OCIijVJvgeec+6glAhGpSVFpOc/OXM9JB3WnT2bboMMREdmPcmT8KC6r4M2FWzjpoB6kp+lgoojErloLPDP71Dl3uJnl440Itu8hwDnnNEa9RNwr8zaxp7icqyZlBx2KiMg+ypHx5/0l28gvKefs0b2CDkVEpElqLfCcc4f7/9NbLhyR/T31+ToO6pXB2AM6BR2KiMg+ypHx58U5G+mWkcrEAZ2DDkVEpEl09rBErRXb8lm6NZ/zxvbBzIIOR0RE4tSKbflMX7aDiyb0JTFB+UZEYpsKPIlaby7cihmcfJBGIBcRkch5+JM1pCUncPnE7KBDERFpMhV4ErXeWrSF8Qdk0jUjLehQREQkTm3PL+aVeZs4b2wfMtulBB2OiEiThVXgmdkBZnacf7uNmemcA4kY5xz3TV/J0q35nDJCrXciEt2UI2PX9vxifvTveZRVVnL14f2CDkdEpFnUW+CZ2XeBF4EH/Em9gVfDXYGZJZrZPDN73b+faWbvmdkK/79Gz5B9CkrK+cEzc/nL28s49eAeXDihb9AhiYjUqik5UvkxWJ+vyuHUuz5lwcY87jxvJNld2gUdkohIswinBe+HwGHAHgDn3AqgawPW8RNgScj9G4APnHODgA/8+yKs3lHA2ff+j3cWb+W3pwzlnotG60KzIhLtmpIjlR8DUFnpuHfaSi55+AvS05J49YeHcc6Y3kGHJSLSbMIp8Eqcc6VVd8wsif2v+VMrM+sNnAo8HDL5TOAJ//YTwFlhRSpx7b2vtnHmPf9jZ0EJT119CNccOUAjZ4pILGhUjlR+DEZuYSnffmIWf31nGace3JOp1x3O0O66ZKGIxJdar4MX4iMz+y3QxsyOB34A/DfM5f8D+BUQej5CN+fcFgDn3BYza0hroMSZvKJS/vjmEp6fvZERvTpw/6Vj6N2pbdBhiYiEq7E58h8oP7YY5xwvzd3En95cQn5xOX846yAuPaSvDiSKSFwKp8C7AbgaWAh8D3iT/Y841sjMTgO2O+fmmNnkhgZmZtcA1wD07avzsOKNc443Fm5hytTF7Coq49qjBvDT4wapS6aIxJoG50jlx5a1fFs+N76yiJlrcxnTtyO3njWCYT3Vaici8SucAq8N8Khz7iHwTgr3pxXV87zDgDPM7BQgDcgws6eBbWbWwz862QPYXtOTnXMPAg8CjBs3LqwuoRIbNuft5abXFvH+ku2M6NWBJ749geE9OwQdlohIYzQmRyo/toCi0nL++cEKHvlkDe3Tkrj9nBGcP64PCbqQuYjEuXDOwfsAL1lVaQO8X9+TnHO/cc71ds5lAxcCHzrnLgWmAlf4s10BvNagiCVmVVQ6nvx8LSf8/WM+XbmT351yIK/8YJKKOxGJZQ3OkcqPkeWc4+1FWznuzo944KPVnDOmFx/+fDIXTuir4k5EWoVwWvDSnHMFVXeccwVm1pSTpG4Hnjezq4H1wHlNWJbEAOcc05ft4E9vLWH5tgKOGNSF284aQd/OOtdORGJec+ZI5ccmWrx5N7e9sYTPVuUwtHs6d100mnHZmUGHJSLSosIp8ArNbIxzbi6AmY0F9jZkJc656cB0/3YOcGzDwpRYtXDjbv745hI+X51Ddue23HfJGE4+qLtObBeReNGkHKn82Dy27ynmzneX8/ycDXRsk8wtZw7nogl9SU4Mp6OSiEh8CafA+ynwgplt9u/3AC6IWEQSFzbkFnHHu8t4bf5mMtulMOX0YVx8yAGkJCnZikhc+SnKkYEpLqvg4U9Wc9/0VZRVVHL1Yf340TGD6NA2OejQREQCU2+B55ybZWZDgSGAAUudc2URj0xi0u69Zdw7bSWP/28tZvCDyQO4dvIAMtKUbEUk/ihHBsM5x9QFm/nzW0vZvLuYE4d34zcnH0h2l3ZBhyYiErhwWvDAS1zD8Eb7Gm1mOOeejFxYEmvKKyp5dtYG/vbecnYVlXLO6N78/ITB9OzYpv4ni4jENuXIFrRyewE3vbaIz1blMLxnBneeP4qJAzoHHZaISNSot8Azs5uByXjJ603gZOBTQMlLAPjfyp384fWvWLo1nwn9MrnptGEc1EsjY4pI/FOObDl7Syu4Z9oKHvx4NWnJifzhrIO4eEJfEjUypojIfsJpwTsXGAnMc85dZWbdCONC5xL/1uws5LY3lvD+km307tSG+y8Zw0kaQEVEWhflyBbwwZJt3Dx1MRt37eWc0b34zSkHkpWeGnRYIiJRKZwCb69zrtLMys0sA+/Cq/0jHJdEsb2lFfzj/eU8+r81pCQm8KuThvDtw/qRlpwYdGgiIi1NOTKCdhaU8NuXF/LuV9sY2LU9z15zKIf2V3dMEZG6hFPgzTazjsBDwBygAJgZyaAkes3fkMfPnp/P6h2FnDe2N788cQhdM9KCDktEJCjKkREyd/0ufvD0XHYVlfKrk4bwncP7ayRmEZEw1Frgmdlhzrn/Adc750qAf5nZ20CGc+7LFotQokJZRSV3f7iSe6etpFt6Ks985xAOG9gl6LBERAKhHBk5zjn+PXM9U6YupnuHNF7+wSSG99R53SIi4aqrBe8uYCzwOTAGwDm3tgVikiizYls+1z8/n0Wb9nDOmF5MOWO4LnsgIq2dcmQEFJdV8PtXF/HCnI0cNTiLf144io5tU4IOS0QkptRV4JWZ2WNAbzO7q/qDzrkfRy4sCVplpWPR5t28tWgrj3y6hvapSfzr0jGcdFCPoEMTEYkGypHNqLisgmlLt3Pv9JUs2rSHHx8zkJ8cN1gjZIqINEJdBd5pwHHAMXjnFUicKy6r4H8rd/L+ku18sGQb2/NLSDA4cXh3bjnzII1YJiLyNeXIJiotr+R/K3cydcFm3l28lcLSCrLSU3n48nEcN6xb0OGJiMSsWgs859xOM3sB6Omce6IFY5IWtD2/mGlLt/PeV9v5dOUOissqaZeSyFFDsjh2aDeOHtqVzHbqHiMiEko5snEqKh0zVufw3y8389aireQVldGhTTKnj+zJ6SN7cki/TJISNZCKiEhT1DmKpnOuwsxOB/7WQvFIhDnnWLo1nw+WbOP9JduZvyEPgF4d23DBuD4ce2A3DumfSWqSLnkgIlIX5cjwVFY65m3YxX8XbOGNhVvYkV9Cu5REjh/WjTNG9eTwgVkaHVNEpBmFc5mEz8zsHuA5oLBqonNubsSikmZXWel47LO1PPrpGjbl7QVgZJ+O/Pz4wRw3rBtDu6frAuUiIg2nHFmHl+du5M53l7Mpby8pSQkcO7Qrp4/sydFDutImRQcSRUQiIZwCb5L//5aQaQ7vvAOJATkFJfz8hQVMX7aDSQM686NjBnLM0K66fp2ISNMpR9bAOce/PlrNn99eyui+Hfn5CYM5flg30jUCs4hIxNVb4Dnnjm6JQCQyvlidw0+enceuwjL+cOZwLj30ALXUiYg0E+XIb6qsdPzprSU89MkazhjZkzvOG6kumCIiLajeAs/MbqppunPulpqmS3SoqHTcO20l/3h/OQd0bsejV47XhWJFRJqZcuT+yioqueGlhbw0dyNXTDyAm08fToIudSAi0qLC6aJZGHI7DW9o6CWRCUeaw+6iMn7477l8unInZ43qya1nj6B9ajhvtYiINJBypK+i0vGjf8/j7cVbuf64wfz42IHqMSIiEoBwumjeGXrfzO4ApkYsImmSNTsLufrxWWzYVcSfvzWC88f1UYIVEYkQ5civ/eH1r3h78VZ+f9owrj68X9DhiIi0Wo1p1mkL9G/uQKTpvlidw7VPz8GAZ75zKBP6ZQYdkohIa9Mqc+Sjn67h8c/W8p3D+6m4ExEJWDjn4C3EGxEMIBHIYv/RwiQKPD97A797ZSF9M9vy6JXjOaBzu6BDEhGJe8qR8O7irfzhja84cXg3fnvKgUGHIyLS6oXTgndayO1yYJtzrjxC8UgDOef45wcr+Mf7Kzh8YBfuvWQMHdpoGGoRkRbSqnPksq35/OTZ+RzcqwP/uGC0BlQREYkC4YxbnARsdc6tAwYBPzCzjhGNSsLinOPOd5fzj/dX8K0xvXnsqvEq7kREWlarzZF7Syv40X/m0i41iYeuGKcLl4uIRIlwCryXgAozGwg8AvQD/h3RqKRezjn+/PYy7pm2kosm9OGv5x5McqKuMyQi0sJabY689Y2vWL6tgL+dP5Ku6WlBhyMiIr5wKoJKv7vJOcA/nHPXAz0iG5bUxTnHbW8s4V8freLSQ/ty21kj1C1GRCQYrTJHvr1oC8/MWM/3juzPkYOzgg5HRERChHMOXpmZXQRcDpzuT1M/wICUlFfw+1cX8fzsjVw5KZubTx+myyCIiASn1eXIdTmF/OrFLzm4dwd+fsKQoMMREZFqwinwrgKuBW5zzq0xs37A05ENS2qyPb+Ya5+aw9z1efz4mIFcf/xgFXciIsFqVTkyv7iM7zwxm4QE456LxpCSpFMDRESiTTgXOv8K+HHI/TXA7ZEMSr5pwYY8vvfUHHbvLeO+S8Zwyoi47wEkIhL1WlOOrKx0XP/cfFbvLOSpb0+gb+e2QYckIiI1aMyFzqWFvTJvI79+aSFZ7VN56fuTGNYzI+iQRESklfnbe8t5f8l2ppw+jEkDuwQdjoiI1EIFXhRzzvH391dw1wcrOKRfJvddMobO7VODDktERFqZ977axj3TVnLBuD5cMSk76HBERKQOES3wzGwtkA9UAOXOuXFmlgk8B2QDa4HznXO7IhlHLCqvqOTGVxfx7KwNnDe2N388Z4QugyAiEkdiJUdu3V3ML19cwPCeGdxy1nCd+y0iEuXqrRjMbJyZvWJmc83sSzNbaGZfNmAdRzvnRjnnxvn3bwA+cM4NAj7w70uIvaUVXPv0HJ6dtYHrjh7IX3SNOxGRqBTvObKi0vHT5+ZRWl7J3ReNJjVJFzMXEYl24bTgPQP8ElgIVDbDOs8EJvu3nwCmA79uhuXGhV2FpVz9xCzmbcjjD2cO57KJ2UGHJCIitYvrHHn/9JV8sTqXv557MP2z2gcVhoiINEA4Bd4O59zURi7fAe+amQMecM49CHRzzm0BcM5tMbOujVx23NlVWMpFD33B6p2F3HfxGE7WSJkiItEubnPkok27+fv7Kzh9ZE/OHds7qDBERKSBwinwbjazh/G6ipRUTXTOvRzGcw9zzm32E9R7ZrY03MDM7BrgGoC+ffuG+7SYtbuojEsfmcHqnYU8csU4jhiUFXRIIiJSvxbPkS2RH8srKvnNywvp1DaFW888SOfdiYjEkHAvdD4USObr7icOqDd5Oec2+/+3m9krwARgm5n18I9M9gC21/LcB4EHAcaNG+fCiDNm7Sku4/JHZ7BiWwEPXD5WxZ2ISOxo8RzZEvnx8c/WsnDTbu65eDQd2iZHYhUiIhIh4RR4I51zIxq6YDNrByQ45/L92ycAtwBTgSvwLgR7BfBaQ5cdT/KLy7jqsVks3ryH+y8dy9FD1GNVRCSGxF2O3JBbxJ3vLufYoV05VacKiIjEnHAKvC/MbJhz7qsGLrsb8IrfrSMJ+Ldz7m0zmwU8b2ZXA+uB8xq43Lgxd/0urn9uPht37eWei0Zz/LBuQYckIiINE1c5sqLS8YsXFpBg8Iez1DVTRCQWhVPgHQ5cYWZr8M4vMMA55w6u60nOudXAyBqm5wDHNiLWuFFeUcndH67knmkr6Z6Rxn++eygT+mUGHZaIiDRcXOXI+6atZMaaXO44byQ9O7YJKgwREWmCcAq8kyIeRSuydmchP31uPvM35HHO6F5MOXM4GWk6v0FEJEbFTY6cvTaXf3ywgjNH9eRbY3oFHY6IiDRSvQWec26dmY0EjvAnfeKcWxDZsOKPc47nZm3glte/IinBuPui0Zw+smfQYYmISBPES47cvbeMnzw7n54d07hVXTNFRGJaQn0zmNlP8C7k2tX/e9rMfhTpwOLJ6h0FXPzQDG54eSEje3fkneuPVHEnIhIH4iVHPva/NWzK28tdF44mXb1KRERiWjhdNK8GDnHOFQKY2Z+Bz4G7IxlYPCgtr+TBj1dx14crSU1K4I9nj+DC8X1ISNCRURGROBHzObK8opJnZ27gyMFZjO7bKehwRESkicIp8AyoCLlf4U+TOsxZt4vfvPwly7cVcOqIHtx8+jC6ZqQFHZaIiDSvmM+RHyzdztY9xdxy5vCgQxERkWYQToH3KDDDvwgrwFnAIxGLKMYVl1XwxzeX8NQX6+iRkcbDl4/jOF3+QEQkXsV8jnz6i3X06JDGMUN1HVYRkXhQZ4FnZgnADOAjvKGgDbjKOTevBWKLOTvyS7jmqdnMW5/HlZOy+cWJQ2ifGk4NLSIisSZecuSCDXmcenBPkhLrPS1fRERiQJ3Vh3Ou0szudM5NBOa2UEwxadnWfL79+CxyCku4/5IxnDyiR9AhiYhIBMVLjuzcPpWCkvKgwxARkWYSzuG6d83sW6Yxk2s1fdl2vnX/Z5RWVPL89yaquBMRaT1iPkd2aZ/CjvzioMMQEZFmEk7/wZ8B7YByMyvG64LinHMZEY0sRjwzYx2/f3URQ7pn8MgV4+jZsU3QIYmISMuJ+RzZpX0qy7blBx2GiIg0k1pb8MzsMP9mlnMuwTmX4pzLcM6lx1LiiqSHPl7N715ZxFGDs3jx2okq7kREWol4ypGj+3Zk9Y5CnpmxLuhQRESkGdTVRfMu//9nLRFIrLn7gxXc9uYSTh3RgwcvH0c7DaYiItKaxE2OvPrw/hw9JIubX1vMjNU5QYcjIiJNVFdVUmZmjwG9zeyu6g86534cubCil3OOO95dxr3TVnHO6F785dyDNfKYiEjrEzc5MjHB+OdFoznr3v/x/WfmMvW6w+jdqW3QYYmISCPVVZmcBrwD7AXm1PDXKt357nLunbaKiyb04Y7zRqq4ExFpneIqR2akJfPw5eMoq6jkO0/MZrsGXRERiVm1tuA553YCz5rZEufcghaMKWrNXpvLvdNXct7Y3vzx7BHE8KBpIiLSBPGYI/tntee+S8ZwzZNzOO2uT7n/0jGMPSAz6LBERKSB6m1+ipfE1VQl5RXc8PJCenZow5Qzhqu4ExGRuMuRRwzK4pUfTqJNSiIXPvgFT36+Fudc0GGJiEgDqH9hmO6dtoqV2wu49eyDNKCKiIjEraHdM5h63eEcOSiLm15bzM9fWEBxWUXQYYmISJhU4IVhzc5C7p++krNG9eToIV2DDkdERCSiOrRJ5qHLx3H9cYN5Zd4mzrnvM9blFAYdloiIhKHepigzSwW+BWSHzu+cuyVyYUWXF+dsoNLBb089MOhQREQkisRzjkxIMH5y3CAO7t2Bnz43n9Pu+pS/nncwJx3UI+jQRESkDuG04L0GnAmUA4Uhf62Cc47Xv9zCpAGd6ZqeFnQ4IiISXeI+Rx49tCtv/Phw+ndtz7VPz+WW/35FaXll0GGJiEgtwjmZrLdz7qSIRxKlFm/ew7qcIn4weUDQoYiISPRpFTmyd6e2vPC9ifzxzSU8+r81zNuwi3suHkOvjm2CDk1ERKoJpwXvMzMbEfFIotR/v9xMUoJx4vDuQYciIiLRp9XkyJSkBKacMZx7Lx7Dim0FnHrXJ0xbtj3osEREpJpaW/DMbCHg/HmuMrPVQAlggHPOHdwyIQanstIxdf5mjhycRce2KUGHIyIiUaI158hTD+7BsJ4ZfP/pOVz12Cx+ePQArj9uMEmJGrdNRCQa1NVF87QWiyJKzViTy5bdxfzmFA2uIiIi+2nVObJfl3a8+sPDmDJ1MfdOW8Wcdbu468LRdM3QueoiIkGr9XCbc26dc24dcGvV7dBpLRdicF6bv4l2KYkcf2C3oEMREZEoohwJacmJ3P6tg7njvJHM35DHKXd9ymerdgYdlohIqxdOf4rhoXfMLBEYG5lwokdZRSVvLNzCicO70yYlMehwREQkOrXKHBnq3LG9ee2Hh5PRJolLH57Bb15eyLY9xUGHJSLSatVa4JnZb8wsHzjYzPaYWb5/fzvesNBxbfHmPeQXl3OsWu9ERKSa1p4jqxvSPZ2p1x3O5ROzeXHOBo766zT+8vZSdu8tCzo0EZFWp64umn9yzqUDf3XOZTjn0v2/zs6537RgjIGYuSYHgPH9OgUciYiIRJvWniNr0j41iSlnDOeDn03mpOHduW/6Ko78yzQe+ng1xWUVQYcnItJqhNNF87dmdo6Z/c3M7jSzsyIdVDSYuSaX/l3a6eLmIiJSl1aZI+vSt3Nb/nHhaN748eGM6tOR295cwjF3TOeF2RuoqHRBhyciEvfCKfDuBa4FFgKLgGvN7N6IRhUw5xyz1+1ifHZm0KGIiEh0a3U5MlzDe3bgiW9P4N/fPYSs9FR++eKXnH3f/1iyZU/QoYmIxLVwCryjgBOdc4855x4DTgEmh7NwM+toZi+a2VIzW2JmE80s08zeM7MV/v+o6wO5o6CEvKIyhvXMCDoUERGJbo3KkbGaHxtj0oAuvPrDw/jnhaPYnLeX0+/+lL+9u4yScnXbFBGJhHAKvGVA35D7fYAvw1z+P4G3nXNDgZHAEuAG4APn3CDgA/9+VFm1vRCAAVntA45ERESiXGNzZEzmx8YyM84c1Yv3rj+KM0b25K4PV3LqXZ8yZ92uoEMTEYk74RR4nYElZjbdzKYDXwFZZjbVzKbW9iQzywCOBB4BcM6VOufygDOBJ/zZngDOanT0EbJqRwEAA7q2CzgSERGJcg3OkbGcH5uqU7sU/nbBKB6/ajx7Sys491+fMWXqYgpLyoMOTUQkbiSFMc9NjVx2f2AH8JiZjQTmAD8BujnntgA457aYWddGLj9i1uUUkpacQPcMDbAiIiJ1akyOjNn82FwmD+nKO9cfyV/fXsrjn61l5ppcnvj2BLLSU4MOTUQk5tXbguec+whYCyT7t2cCc51zH/n3a5MEjAHud86NBgppQHcTM7vGzGab2ewdO3aE+7RmkV9cToc2yZhZi65XRERiSyNzZMzmx+bUPjWJ/zvzIB67ajxrdhZy7r8+Y31OUdBhiYjEvHoLPDP7LvAi8IA/qTfwahjL3ghsdM7N8O+/iJfQtplZD3/ZPfAuCvsNzrkHnXPjnHPjsrKywlhd8yksraBdSjiNmyIi0po1MkfGbH6MhKOHdOXp7xxCXlEZ3/rXZxplU0SkicI5B++HwGHAHgDn3Aqg3m4jzrmtwAYzG+JPOhbv3ISpwBX+tCuA1xoYc8QVlZTTJiUx6DBERCT6NThHxnJ+jJSxB3TihWsnkmjG+Q98zsw1uUGHJCISs8Ip8Eqcc6VVd8wsCQj3SqU/Ap4xsy+BUcAfgduB481sBXC8fz+q7Ckuo32qWvBERKRejc2RMZkfI2lwt3Re/P5EstqncvFDX/DP91dQVlEZdFgiIjEnnCrmIzP7LdDGzI4HfgD8N5yFO+fmA+NqeOjYsCMMwMZde5k0oEvQYYiISPRrVI6M1fwYab07teWVHxzGTVMX8ff3l/Ph0m387YJRumyRiEgDhNOCdwPeaF8Lge8BbwI3RjKoIBWXVbB1TzF9M9sGHYqIiES/VpUjW0KHtsn888LR3HPxaNblFnHKPz/h8f+tobIy3M5DIiKtW70teM65SjN7FXjVORe7w3WFaVPeXpyDvp3bBB2KiIhEudaWI1vSaQf3ZEJ2Jr9+6Uum/Pcr3luyjdvPOZg+OgArIlKnWlvwzDPFzHYCS4FlZrbDzBp7XbyYsGKbd5HzAzrrIuciIlKz1pojW1rXjDQevXI8fzpnBPPW53Hs3z7iz28vJb+4LOjQRESiVl1dNH+KNzLYeOdcZ+dcJnAIcJiZXd8SwQVh7vpdpCQmMKxHRtChiIhI9PoprTBHBsHMuGhCX97/2VGcOqIH909fxdF3TOffM9ZToW6bIiLfUFeBdzlwkXNuTdUE59xq4FL/sbg0e20uI3p3IC1Zl0kQEZFatcocGaSeHdvw9wtG8doPD6Nfl3b89pWFnHrXJ3yyQj1jRURC1VXgJTvndlaf6J9jkBy5kIJTXFbBok17GHdAp6BDERGR6NbqcmS0GNmnI89/byL3XTKGwtJyLntkJt9+fBYbcouCDk1EJCrUVeCVNvKxmLViWwGlFZWM7NMx6FBERCS6tbocGU3MjFNG9OC964/iNycPZeaaXE67+1M+WLIt6NBERAJXV4E30sz21PCXD4xoqQBb0uqd3gArut6OiIjUo9XlyGiUlpzI944awBs/Ppzendpw9ROz+dNbSyjXBdJFpBWr9TIJzrlWdxLa2p1e944DOmsIZhERqV1rzJHR7IDO7Xjp+5P4w+tf8cBHq5m3Lo+7LhpN9w5pQYcmItLiwrnQeauxZmcBvTq20QArIiIiMSYtOZHbzh7BPy8cxaLNuzn1rk/4eLkGYBGR1kcFXoi1OUVqvRMREYlhZ47qxdTrDqNz+xQuf3Qmf357KWXqsikirYgKvBAbclXgiYiIxLqBXdN57YeHc9GEvtw/fRXnP/C5RtkUkVZDBZ6vsKScnMJS+mSqwBMREYl1bVIS+dM5I7jn4tGs3FbAKXd9wpsLtwQdlohIxKnA823Y5R3Z66sCT0REJG6cdnBP3vzJEfTPas8PnpnLbW98FXRIIiIRpQLPt2V3MQA9NOKWiIhIXOmT2ZYXvjeRSw7py0OfrOEtteSJSBxTgefLLfCuS9u5XWrAkYiIiEhzS0lKYMoZwxnRqwO/e3URO/JLgg5JRCQiVOD5cgu9Ai+zfUrAkYiIiEgkJCcmcOf5IykoKed3ryzEORd0SCIizU4Fnm9nYQnJiUZ6aq3XfhcREZEYN7hbOr84YTDvfrWNF+dsDDocEZFmpwLPN3NNLkO6p2NmQYciIiIiEXT14f05tH8mv39tEUu37gk6HBGRZqUCD9i2p5h56/M4cVj3oEMRERGRCEtMMO66aDQZacl8/+m57CkuCzokEZFmowIPePerbQCceJAKPBERkdaga3oa91w8hvW5RfzqhS91Pp6IxA0VeMB7X22jX5d2DOraPuhQREREpIVM6JfJDScN5e3FW3nok9VBhyMi0ixafYFXWemYt24XkwZ01vl3IiIircx3jujHScO78+e3lzFjdU7Q4YiINFmrL/DW5hSSX1LOwb07BB2KiIiItDAz46/nHUzfzLZc9595bM8vDjokEZEmafUF3sJNuwEY0atjsIGIiIhIINLTkrn/0jHkF5dx3b/nUV5RGXRIIiKN1uoLvEWbdpOSlMCgbjr/TkREpLUa2j2DP50zgplrcrnrgxVBhyMi0mitvsBbtq2AQV3bk5zY6jeFiIhIq3b26N6cM7oX901fxbKt+UGHIyLSKK2+qlm1vYCBGj1TREREgBtPG0ZGm2RuePlLKip16QQRiT2tusArLClnU95eBmapwBMRERHIbJfC7087kHnr83j6i3VBhyMi0mCtusBbunUPgM6/ExERkX3OGtWLIwZ14a/vLKOwpDzocEREGiRiBZ6ZDTGz+SF/e8zsp2aWaWbvmdkK/3+nSMVQn/8u2EJKUgITB3QJKgQREWmFYiFHtmZmxvcnD6CgpJyPl+8IOhwRkQaJWIHnnFvmnBvlnBsFjAWKgFeAG4APnHODgA/8+y2utLySqQs2c/yB3ejQJjmIEEREpJWK9hwpMCE7k05tk3l78dagQxERaZCW6qJ5LLDKObcOOBN4wp/+BHBWC8Wwn4+W7yC3sJRzxvQKYvUiIiJVoi5HCiQlJnDcgd34cMl2Sst1XTwRiR0tVeBdCPzHv93NObcFwP/ftYVi2M/rX24ms10KRw7OCmL1IiIiVaIuR4rnxOHdyS8pZ9ba3KBDEREJW8QLPDNLAc4AXmjg864xs9lmNnvHjubv/z577S4m9u+s69+JiEhgGpMjI50f5WsT+mdi5v1mEBGJFS1R3ZwMzHXObfPvbzOzHgD+/+01Pck596BzbpxzblxWVvO2sm3ZvZdNeXsZe4DOXRcRkUA1OEdGMj/K/jLSkhnSLZ0561XgiUjsaIkC7yK+7noCMBW4wr99BfBaC8Swn6ojceOyVeCJiEigoi5Hyv7GHNCJeet2UamLnotIjIhogWdmbYHjgZdDJt8OHG9mK/zHbo9kDDWZvyGP1KQEDuyR0dKrFhERAaI3R8r+RvXpSH5JOetyi4IORUQkLEmRXLhzrgjoXG1aDt6IYYFZtjWfId3Tdf6diIgEJlpzpOyvZ4c2AGzfU0y/Lu0CjkZEpH6tssJZujWfwd3Sgw5DREREolxWeioAOwpKAo5ERCQ8ra7AyykoYWdBCUO7q8ATERGRunVpnwLAznwVeCISG1pdgbdw024AnX8nIiIi9erUNoXkRGObCjwRiRGtrsD7dMVOUpISdIkEERERqVdCgtGrYxvW52iQFRGJDa2uwPtkxU7GZ3ciLTkx6FBEREQkBvTt3I51uYVBhyEiEpZWVeBt31PMsm35HD5QF4YVERGR8ByQ2ZZ1OUU4p2vhiUj0a1UF3rwNeQAc2j8z2EBEREQkZvTs2Ib84nL2llUEHYqISL1aVYG3ZqfXvWJA1/YBRyIiIiKxon2qd1pHYYkKPBGJfq2rwNtRSJf2KWSkJQcdioiIiMSINilJAOwtVYEnItGvdRV4Owvp16Vd0GGIiIhIDGmX4rfglZYHHImISP1aVYG3akeBCjwRERFpkLapXgtekQo8EYkBrabA21lQQk5hKYO7pQcdioiIiMSQtn4LXpG6aIpIDGg1Bd7ybfkAKvBERESkQaoKPA2yIiKxoNUUeCu2FQAwpLsKPBEREQlf2xR10RSR2NFqCrx563fRpX0KXdNTgw5FREREYsjXg6yoBU9Eol+rKPCcc3y2KoeJA7pgZkGHIyIiIjGkjV/g7VULnojEgFZR4K3eWcj2/BImDegcdCgiIiISY77uoqkWPBGJfq2iwPtsVQ4AE/urwBMREZGGSUww0pITVOCJSExoFQVem+REjhycxQGd2wYdioiIiMSgowZn0SdTvyNEJPolBR1ASzh3bG/OHds76DBEREQkRj1w2bigQxARCUuraMETERERERFpDVTgiYiIiIiIxAkVeCIiIiIiInFCBZ6IiIiIiEicUIEnIiIiIiISJ1TgiYiIiIiIxAkVeCIiIiIiInFCBZ6IiIiIiEicUIEnIiIiIiISJ1TgiYiIiIiIxAkVeCIiIiIiInFCBZ6IiIiIiEicUIEnIiIiIiISJ8w5F3QM9TKzHUAhsDPoWBqoC4o50mItXlDMLSHW4gXFXOUA51xWMy8zbvn5cV0TFhFr+53ijbxYi1nxRlasxQuxF3ND4g0rR8ZEgQdgZrOdc+OCjqMhFHPkxVq8oJhbQqzFC4pZghFr76HijbxYi1nxRlasxQuxF3Mk4lUXTRERERERkTihAk9ERERERCROxFKB92DQATSCYo68WIsXFHNLiLV4QTFLMGLtPVS8kRdrMSveyIq1eCH2Ym72eGPmHDwRERERERGpWyy14ImIiIiIiEgdVOCJiIiIiIjEiagq8Mysj5lNM7MlZrbYzH5SwzyTzWy3mc33/24KItZqMa01s4V+PLNreNzM7C4zW2lmX5rZmCDi9GMZErLt5pvZHjP7abV5At/GZvaomW03s0Uh0zLN7D0zW+H/71TLc08ys2X+9r4h4Jj/amZL/ff9FTPrWMtz69yHWjjmKWa2KeT9P6WW57b4dq4l3udCYl1rZvNreW5Q27jG77Vo3Z/riDeq9+XWzszSzGymmS3w37f/q/b4L8zMmVmXWp5f4/sW7n7akvHWlcfC/f5qzpib+p3Z0ts4nHhr+x5oyOttyXj9+Vp0H25qzEHsx3V97szsR/7+udjM/lLL86NiHw4n3mjah8OJ15+n+fZh51zU/AE9gDH+7XRgOTCs2jyTgdeDjrVaTGuBLnU8fgrwFmDAocCMoGP240oEtuJdNDGqtjFwJDAGWBQy7S/ADf7tG4A/1/KaVgH9gRRgQfV9qIVjPgFI8m//uaaYw9mHWjjmKcAvwth3Wnw71xRvtcfvBG6Ksm1c4/datO7PdcQb1ftya//z80t7/3YyMAM41L/fB3gH74LoNb43tb1v4eynQcQbspz98lg431/NHXNTvzNbehuHGW+tv8citY2bEm8Q+3BzxNzS+3Ed8R4NvA+k+o91jfJ9OJx4o2kfrjfe5t6Ho6oFzzm3xTk317+dDywBegUbVbM4E3jSeb4AOppZj6CDAo4FVjnn1gUdSHXOuY+B3GqTzwSe8G8/AZxVw1MnACudc6udc6XAs/7zIq6mmJ1z7zrnyv27XwC9WyKWcNWyncMRyHauK14zM+B84D+RjqMh6vhei8r9ubZ4o31fbu38/FLg3032/6pGUfs78KuQ+w0Rzn7aYM0Yb4vlsXpirk9dn+UgtnF9z23x32NN3L51icj2hWaNuUX24zri/T5wu3OuxJ9vew1Pj6Z9uN54o2wfDmf71qXB2zeqCrxQZpYNjMarfqub6Dd/vmVmw1s2sho54F0zm2Nm19TweC9gQ8j9jURH4Xohtf8YjrZtDNDNObcFvA8u0LWGeaJ1WwN8G68ltyb17UMt7TrzuuI9WktXgGjczkcA25xzK2p5PPBtXO17Ler35zq+h2NpX241zCzRvC7K24H3nHMzzOwMYJNzbkE9T6/tfQtnPw0i3io15bH6vr8araaYw1xnXZ/lFt3GYcYbuoxsvvk9EJFt3MR4W3wfboaYq7TYflxLvIOBI8xshpl9ZGbja3hqNO3D4cQbuoxsgt2Hw4232fbhqCzwzKw98BLwU+fcnmoPz8Vrwh4J3A282sLh1eQw59wY4GTgh2Z2ZLXHrYbnBHp9CjNLAc4AXqjh4WjcxuGKum0NYGa/A8qBZ2qZpb59qCXdDwwARgFb8Lo9VheN2/ki6m69C3Qb1/O9VuvTapjWItu5tnhjbF9uVZxzFc65UXitqxPM7GDgd0A451G3+PvWxHhry2PhfH81Wg0xHxTmOgP5LDchXqDW74GIbeMmxhvId08zbOMW3Y9riTcJ6ITXnfCXwPN+r5j9Qq1pcc0VV22aEC8QNftwuPE22z4cdQWemSXjvRHPOOderv64c25PVfOnc+5NINlqOWm8pTjnNvv/twOv4DVjh9qId05Bld7A5paJrlYnA3Odc9uqPxCN29i3raprq/+/pibuqNvWZnYFcBpwiXOuxi/DMPahFuOc2+Z/QVUCD9USS1RtZzNLAs4BnqttniC3cS3fa1G7P9f2PRxr+3Jr5ZzLA6bjdevpBywws7V4+89cM+tew3Nqe9/C2U9bPF7fN/JYmN9fzRnzSc3wndmS2zjceGv9HmiJbdyYeIPchxsbsy+Q/Tg0Xrz982W/i+FMoBKo/rsvavbhMOONmn043Hibcx+OqgLPr2YfAZY45/5Wyzzdq6peM5uA9xpyWi7Kb8TTzszSq27jDUSwqNpsU4HLzXMosLuqqTVAtbZ2RNs2DjEVuMK/fQXwWg3zzAIGmVk//6jYhf7zAmFmJwG/Bs5wzhXVMk84+1CLsf3PDz27lliiajsDxwFLnXMba3owyG1cx/daVO7PtcUbi/tya2JmWeaPbGpmbfA+E/Occ12dc9nOuWy8HxljnHNbqz23rvctnP20ReMN8Y08Fub3V3PGvLQZvjNbchuHFW9dv8citY2bGG+L78NNjTlEi+3HtcWL11PrGH/6YLxBVHZWe3rU7MPhxBtN+3CY8TbvPuyaeQSZpvwBh+M1934JzPf/TgGuBa7157kOWIw3es8XwKSAY+7vx7LAj+t3/vTQmA24F2/0oYXAuIBjbotXsHUImRZV2xjvy24LUIaX5K8GOgMfACv8/5n+vD2BN0OeewreaEmrqt6PAGNeiddnvWp//lf1mGvbhwKM+Sl/P/0S70ulR7Rs55ri9ac/XrX/hswbLdu4tu+1qNyf64g3qvfl1v4HHAzM89+3RdQwmiwhI7SF+77Vtp8GGa9//xt5zJ9e4/dXJGOubZ3hfpZbehuHE29t3wOR3MZNjLfF9+Fm2idadD+uI94U4Gl/2lzgmCjfh+uNN8r24XDibdZ92PwnioiIiIiISIyLqi6aIiIiIiIi0ngq8EREREREROKECjwREREREZE4oQJPREREREQkTqjAExERERERiRMq8ERERERikJlVmNl8M1tkZi+YWdtmWOZkM5sUcv9aM7u8kcsqqGX678xssZl96cd/SD3LedzMzvVvP2xmw+qZf1/MZnalmfVsYNz1rqPa/FPMbJOZ3RLm/A+Y2WFmNt3MxjVgPdeb2Xozuyfc50jrlBR0ACIiIiLSKHudc6MAzOwZvGva7ruos5klOucqGrjMyUAB8BmAc+5fzRLp1zFNBE7Du4h9iZl1wbtOWFicc98JY57QmK/Eu/7Y5uZcRw3+7py7I8x5DwF+0NAVOOf+bma7gLCLQmmd1IInUaehR/bCWF5HM2vwF2kQWiJWM/usAfNmm9leM5tfy+NTzOwXzRDTeWY2w3/PV5rZzWbWxn//S/0fACIiUrtPgIF+C9w0M/s33sWcMbNXzWyOn1uvqXqCmZ1kZnPNbIGZfWBm2XhF4vX+9+8Rod/zZjbQzN73559rZgPMrL3/3LlmttDMzqwnzh7ATudcCYBzbqdzbrO//JvMbJbfIvmgmVn1J4e2eplZgZnd5sfzhZl186dPMbNf+K1+44Bn/Ndzqpm9ErKs483s5casoy5mlmVm7/nb5AEzW1eVx8zsQGB5aOFtZglm9oSZ3erfv9rMlvtxPKQWO2koFXgSVaod2TsYOA7Y0MTFdqQRR8oixTy1ffY60sBY61neNzjnJtU/135WVR0hjgQzuwL4NfAt/z0fBRQ556qOTId91FVEpDUysyTgZPyCDpgA/M45V9XN8NvOubF4xc6PzayzmWUBD+F9944EznPOrQX+hdcaNco590m1VT0D3OvPPwnYAhQDZzvnxgBHA3fWVJiFeBfo4xcw95nZUSGP3eOcG++cOwhog/d7oC7tgC/8eD4Gvhv6oHPuRWA2cImfT94EDvRfO8BVwGNNWUctbgY+9LfJK0DfkMdOBt4OuZ+Et12XO+duNK876e+BQ4HjgaFhrE9kPyrwJNrUdWTvUjOb6R+Fe8DMEv3pNR6ZDHE7MMB/3l9rW5bfWrXUvL73i8zsGTM7zsz+Z2YrzGyC/9yq+Z4wr8XpRfPPe6gjxmwzW2Jm9wFz8ZJbTXHvF6v/vEVVL8Q/IjmlluXVuO7q/KORVc9/yF//u2bWJpw3yLwW1mVm9j4wJGR6ba/99/72es/M/mMhLX5mloHXneh859xG/z0vcM79NZxYRERauTbm9bCYDawHHvGnz3TOrQmZ78dmtgD4AugDDMIrID6ums85l1vXiswsHejlnHvFn7/YOVcEGPBHM/sSeB/oBdTayuWcKwDGAtcAO4DnzOxK/+GjzevNsRA4Bhhez+svBV73b88Bsuua2TnngKeAS82sIzAReKs51+E7HHjWX+fbwK6Qx05k/wLvAWCRc+42//4E4CPnXK5zrgx4IYz1iexHBZ5EmxqP7PldGi4ADvOPwlUAl/jP+caRyWrLvAG/Fco598t6ljUQ+CdwMN5Rs4vxvqh/Afw2ZJlDgAf9Fqc9wA/qWW7Vc550zo12zq2rJe79Yq1nW+1bHtC2nnXXZBDekdjhQB7wrXrmx8zGAhcCo4FzgPH+9Bpfu3ldXL4VMn/18wbOBmY451bXt24REfmGvX6+GOWc+5FzrtSfXlg1g5lNxusNM9FvhZoHpOEVZq4B66qtVe4SIAsY63//b/OXXyvnXIVzbrpz7mbgOuBbZpYG3Aec65wbgde6WOdygDK/aAMv74QztsRjwKXARcALzrnyCKyjxm3lHwzuWHXg2vcZXmFb9Vrrav0UCYsGWZGo4pwr8IuII/C6ejxnZjcA7fGO+M3ye360Abb7T/uxmZ3t3646MplTx2qOrWVZHwNrnHNV5ywsBj5wzjn/aGJ2yDI2OOf+599+GvgxXjeV2mIEWOec+yLkfk1xb60j7upCl1fba6rLGufcfP92uEcljwBe8Y/aYmZT61l/JvCac26vP/9/qy1vODAfERGJlA7ALudckZkNxWu5A/gcuNfM+jnn1phZpt+Klw9kVF+Ic26PmW00s7Occ6+aWSqQ6C9/u3OuzMyOBg6oKxgzGwJUOudW+JNGAev4upjbaWbtgXOBF5vywn35QHrI69hsZpuBG/G6QEbCp8D5wJ/N7ASgkz/9aGBatXkfAY4EXvB/E8wE/m5mnfzYv8XXXW9FwqICT6KOf+LxdGC6X1hdgdey94Rz7jeh81Y7MllkZtOp/4if1bKsbKAkZFJlyP1K9v+8VD/q6WpbbojajqjWFXc5+7e0h85TGHK7vnXXJPS1VuAVZeGo6Yhvbdv0+nqWVdiA9YqISMO9DVzrd6FchtdNE+fcDv/0gJfNO497O17B81/gRfMGS/lRtWVdBjxg3uUAyoDz8M4f+6+ZzcY7YLe0nnjaA3f7XSTLgZXANc65PDN7CK+YWQvMasqLDvE48C8z24uXc/f6MWc5575qpnVU93/Af8zsAuAjvHMV8/HOv/tG0eqc+5uZdcDrPnoJ8EdgBt456F8BuyMUp8QpFXgSVeo4svcB8JqZ/d05t93MMvGOyNV2ZDLUfkfv6lhWQ/Q1s4nOuc/xunl8Wtty/e6Y1dUWd/VYtwFd/e6bBXgnnL/NNzVk3U3xMfC4md2O9/1xOt75A7Vt00/xfgz8yZ//VLxuN1XeBJ71n7fNPyJ8uXMudB4REamBc659DdOm4x0krbpfgldY1PT8t6h2DppzbjneaQpVPgl5bAXeuXHVTWxAfHPwBmipaf4b8VrWqk+/MuT25JqW7w+o8qJ/e0rI9JeAl6ot8nD2z0XV11fvOuqxGzjROVdu3uBxRzvvkhCTgH0HPqut5+aq22b2b+fcg+YNnvMK3kFukbCpwJNoU9uRvZ1mdiPwrn+ksQz4IbUcmQzlnMsxb6CURcBb/nl4NS2rId0jlwBXmNkDwArgfr9Yq2m5NRVZtR1RrSnWW/CO5K2hliOjzrmvGrDuRnPOzTWz5/CO0q7DT/y1rd8594XfjXOBP/9sQo5EOudmmdkU4B1/UJYkvC6vIiIizc7M5uD1Hvl5My62ALjGzDKcczfhjZr5vJ8PS/FH3nTeqJrhmGJmx+H12nkXeNWP/Xq8y1hUL1hF9mNfnzcqIuHwu3K+7rxhnONac7xWM2vvn1vZFq8F8Brn3NwGPH8tMM45t7OxMYiIiIi0FhpFU0TqUgF0sFoudB6mB/3nzwVeCre4M/9C50Ay3jmQIiIiIlIPteCJiIiIiIjECbXgiYiIiIiIxAkVeCIiIiIiInFCBZ6IiIiIiEicUIEnIiIiIiISJ1TgiYiIiIiIxAkVeCIiIiIiInFCBZ6IiIiIiEic+H/1Noo43mJIQgAAAABJRU5ErkJggg==\n", "text/plain": ["
"]}, "metadata": {"needs_background": "light"}, "output_type": "display_data"}], "source": ["fig, ax = plt.subplots(ncols=2, figsize=(15,6))\n", "\n", "fig.suptitle('Temperature and salinity at depth with increasing pressure', fontsize=16)\n", "\n", "ocean_data.temp.plot(y='n_levels', ax=ax[0], yincrease=False, label='T')\n", "ocean_data.psal.plot(y='n_levels', ax=ax[1], yincrease=False, label='Salinity')\n", "\n", "ax[0].legend()\n", "ax[1].legend()\n", "\n", "ax[0].set_xlabel('Sea temperature in [$\\deg C $]')\n", "ax[0].set_ylabel('Depth from surface in number of levels')\n", "ax[1].set_xlabel('Practical Salinity in [g/kg]')\n", "ax[1].set_ylabel('Depth from surface in number of levels')"]}, {"cell_type": "markdown", "id": "c9a95e9d", "metadata": {}, "source": ["From the above plots, it can be seen that temperature and salinity look homogeneous in the mixed layer.\n", "Below the mixed layer is the thermocline which is a region where the temperature gradient is very steep. The termocline is between $300$ and $1000\\,m$ and is steepest for example in the Artic where fresh cold water overlies warmer watter."]}, {"cell_type": "markdown", "id": "fcb5d6f3", "metadata": {}, "source": ["### Temperature-salinity relationship\n", "Next, a temperature-salinity profile can be plotted that is independent of depth."]}, {"cell_type": "code", "execution_count": 74, "id": "b322de7a", "metadata": {}, "outputs": [{"data": {"text/plain": ["[]"]}, "execution_count": 74, "metadata": {}, "output_type": "execute_result"}, {"data": {"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZkAAAEKCAYAAADAVygjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAA73ElEQVR4nO3dd3xUZfb48c9JIySEmtAJCVWK9N5BEUXXXrCjInZl1bXs7nfXLf5cy7oo2BBQdBF3LVgRC6AUBaQXEQgdQkdCAoEkcH5/3Bsds5OZSZjJTJLzfr3mlZl773PvyXDN8XnuU0RVMcYYY0IhKtwBGGOMqbgsyRhjjAkZSzLGGGNCxpKMMcaYkLEkY4wxJmQsyRhjjAmZcpFkRKSJiMwRkXUislZE7vNyzCARyRKRFe7rT+GI1RhjzC9iwh1AgAqAB1R1mYgkAUtF5EtV/aHIcfNU9YIwxGeMMcaLclGTUdXdqrrMfZ8NrAMahTcqY4wx/pSXmszPRCQN6Aws8rK7t4isBDKBB1V1ra9zJScna1paWtBjNMaYimzp0qUHVDUlkGPLVZIRkWrAe8AYVT1SZPcyoKmq5ojIcOADoKWXc4wGRgOkpqayZMmS0AZtjDEVjIhsC/TYctFcBiAisTgJZqqqvl90v6oeUdUc9/0MIFZEkr0cN0FVu6lqt5SUgBKxMcaYUioXSUZEBJgErFPVZ4s5pr57HCLSA+d3O1h2URpjjCmqvDSX9QWuB1aLyAp32++BVABVfRm4HLhDRAqAXGCE2hTTxhgTVuUiyajqfED8HDMeGF82ERljjAlEuWguM8YYUz5ZkjHGGBMylmSMMcaEjCUZUyo7fzrG05//SObh3HCHYoyJYJZkTKnMWrePF+Zs4nj+yXCHYoyJYJZkTIkdyDnBy99som2D6qQnJ4Y7HGNMBLMkY0pk0eaDXDhuPgeP5vHkZR1wx78aY4xX5WKcjAmP3Vm5/Lgnm417s9mwN4eN+3JYvfMwTesk8vboXpzZuEa4QzTGRDhLMpXYiYKTHMkt4MjxfLJy8zl8LI8DOXks336YBRkH2H7o2M/HpiRVoVW9atw+sDl3Dm5BtSp26xhj/LO/FBWAqnI8/xRZufkcPHqCgzl5HDqax4GcExw66rz/6Vgeh4/lO6/cPLJy8zmef8rr+ZKqxNCreR1u6ptGu4Y1aFWvGjUT4sr4tzLGVASWZCJAXsGp/0kKWbn5v3plH8/nRMEp8tzX8YKTZB8vIPt4AUdy8yk45X2atpgooVZiHLUSYqmZEEfTOgl0SqhJjYRYalSNpXp8DNWrxlI9PpZaiXHUToijYc14YqLtcZ0x5vRZkikjR08UsH5vNuv3/PLac+Q4B3JOkH28oNhyCXHR1KgaS1J8DPGx0cRFRxEXE0VSfAzNU6qRFB9DUryzv0bVWOokxlGnWhVqJ8aRnFiF6lVj7OG8MSZsLMkEmaqyfm82a3YdYePebDbuy2HD3mx2/vTLoMXEuGha1U+ifaMaTlLwTAzV4qiVGOfWMmKJi7EahTGm/LIkEwSqyupdWcxYvYfP1uxm20HngXlcdBTNUhLpnFqLq7o14YwG1TmjfhKNalYlKspqF8aYis+SzGmas34f//fBGnb+lEt0lNCneR1uG9CcHum1SauTYM82jDGVmiWZUlJVXp23mSc++5HW9ZJ46vIODG1Tj1qJ1gvLGGMKWZIpBVXloXdX8c7SnZx/ZgOevqIDCXH2VRpjTFHl5i+jiGwFsoGTQIGqdiuyX4DngOHAMWCkqi4LRSwb9+XwztKd3Nw3nf+7oI313jLGmGKUmyTjGqyqB4rZdx7Q0n31BF5yfwbdzp+cB/sXdGxgCcYYY3yoSE+lLwLeUMdCoKaINAjFhTIPHwegfvX4UJzeGGMqjPKUZBT4QkSWishoL/sbATs8Pu90twVd4fT23289FIrTG2NMhVGekkxfVe2C0yx2l4gMKLLfW7vV/8y1IiKjRWSJiCzZv39/qQLp3awOqbUTeGvR9lKVN8aYyqLcJBlVzXR/7gOmAz2KHLITaOLxuTGQ6eU8E1S1m6p2S0lJKVUsUVHCiB5NWLTlEJv355TqHMYYUxmUiyQjIokiklT4HjgHWFPksI+AG8TRC8hS1d2hiunyro2Ji4nizx+t5WQxk1MaY0xlVy6SDFAPmC8iK4HFwKeqOlNEbheR291jZgCbgQzgVeDOUAZUNymev1zYjnkbD/CvLzeE8lLGGFNulYsuzKq6GejoZfvLHu8VuKss47q6Ryorth9m/JwMOjSuwTnt6pfl5Y0xJuKVl5pMxPrLRe3o0LgGY/6zwnqbGWNMEZZkTlN8bDQTb+hG/erxjJy82BKNMcZ48Jlk3IfoPUXkUhG5xH1vQ9yLqFs9nrdH96Je9XhutERjjDE/KzbJiMg5wEbgMZz5wM4H/gJsdPcZD4WJpn4NJ9Es3mKJxhhjfNVkngPOVtXzVHWU+zoXGOruM0XUrR7P27f2okGNeEa+tpiFmw+GOyRjjAkrX0kmBmeAY1G7gNjQhFP+1a0ez7TRvWhYsyo3vfY9CzKKm8/TGGMqPl9JZjLwvYg8LCLXuK+HgUXApLIJr3yqmxTPtFt7kVo7gZtf/55vNpRu+hpjjCnvik0yqvoEcC3OnGC9gT7u+2vdfcaHlKQqTBvdi2Yp1bh1yhK+WLsn3CEZY0yZE2cMY+XUrVs3XbJkSUivcfhYHjdOXszqXVn849IOXNm9if9CxhgTwURkadGFI4vjq3dZNRH5q4isFZEsEdkvIgtFZGTQIq0EaibE8datvejbIpmH3lvFC3MyqMyJ3RhTufh6JjMVZy6wYThdl58HrgcGi8j/K4PYKozEKjFMurE7F3ZsyNOfr+cxm1TTGFNJ+Eoyaar6uqruVNVngQtVdSNwE3Bp2YRXccTFRDH2qk7c2j+dKd9t466pyziefzLcYRljTEj5SjJHRaQfgIj8BjgEoKqn8L5AmPEjKkr4w/lt+eP5bZi5dg/XTlzEvuzj4Q7LGGNCxleSuR14VkQOAw8D9wCISArwQuhDq7hG9W/GC9d0YW1mFheOW8DKHYfDHZIxxoSEry7Mq1S1h6rWVNV+qrrB3b5fVZ8vuxArpvM7NOC9O/oQEy1c8cp3vLNkR7hDMsaYoCvRLMwi8kmoAqmM2jWswUd396Nb01r87t1V/PGD1facxhhToZR0qv9GIYnCDxFpLSIrPF5HRGRMkWMGuV2tC4/5UzhiLanaiXG8cXMPbhvQjH8v3M4lL37Lpv054Q7LGGOCoqRJZnlIovBDVderaidV7QR0BY4B070cOq/wOFX9a5kGeRpioqN4dHgbJo/sxp6sXH4zbj7vLfU2bZwxxpQvAScZEYkFxolI3RDGE4izgE2qui3McQTdkDPq8dl9A2jfqAYPvLOS+/+7gqMnCsIdljHGlJqvEf8vi0g7930NYCXwBrBcRK4uo/i8GQFMK2ZfbxFZKSKfFcZe3tSv4Uyued9ZLZm+fBe/GT+fHzKPhDssY4wpFV81mf6qutZ9fxOwQVXPxGmueijkkXkhInHAhcA7XnYvA5qqakdgHPBBMecYLSJLRGTJ/v2ROTtydJTw26GtmDqqJznHC7j4xQW8+d1Wm47GGFPu+EoyeR7vh+L+0VbVcE4nfB6wTFX3Ft2hqkdUNcd9PwOIFZFkL8dNUNVuqtotJSUl9BGfhj7Nk5lxX396N6vD/324ljunLiMrNz/cYRljTMBifOw7LCIX4CxS1he4BUBEYoCq3gqIyEcBXPOQqo4sYZyFrqaYpjIRqQ/sVVUVkR44CbTcL02ZXK0Kr43szqvzNvP05+tZtXMe46/pTOfUWuEOzRhj/PKVZG7DmRSzPjDGowZzFvBpMWXaAKN8nFMo5WwBIpKAU6O6zWPb7QCq+jJwOXCHiBQAucAIrSDtS1FRwm0Dm9M9vTb3vLWcK17+jofObc2ofs2IirIZfowxkavY9WRE5G5VHV+ik4lcqar/Pd1jykpZrCcTbFnH8nn4vVXMXLuHIWfU5ZkrOlI7MS7cYRljKpGgrCcD3FzSC3tLHiISJSLVfR1jAlcjIZaXruvCXy5sx/yNBxj+3DwWbzkU7rCMMcarkg7GDIiIvCUi1UUkEfgBWC8ivwvFtSojEeHGPmm8f2cf4mOjGDHhO57+/EfyCk6FOzRjjPkVX0mmgzt9S9FXtoj4G7jRVlWPABcDM4BUnAXPTBC1b1SDT+7tz+VdG/PCnE1c8uICNuzNDndYxhjzM19JZrWqVvfySlLV6j7KgdN9OBYnyXyoqvlAhXgIH2mqVYnhqcs7MuH6ruzJOs4F4+Yzaf4WTtnKm8aYCBCS5jLgFWArkAjMFZGmgA1bD6Fz2tVn5pgBDGiZzN8++YHrJi1i1+HccIdljKnkfCUZb6PqA6Kqz6tqI1Udro5twODSns8EJiWpCq/e0I1/XHomK3YcZti/5jJt8XabKcAYEza+xsk0EJFiFydT1XuL2+djmv1yMzNyeSUijOiRSt8WyTz07ioefX81M1bv5h+XdaBRTa9jaI0xJmR8JZnTGUBy1ON9PHABsO40zmdKqEntBKaO6snURdt44rMfGfavufzx/DZc1b0JIjaA0xhTNoodjBnUi4hUAT5S1WEhv1gJlMfBmKWx/eAxHnpvJQs3H2JAqxSeuPRMq9UYY0otWIMxgykBaFZG1zJFpNZJ4K1RvfjrRe1YsvUQQ5/9honzNlNw0sbVGGNCK1SDMVeLyCr3tRZYDzwXimuZwERFCTf0TuPzMQPomV6bv3+6jgvHL2DFjsPhDs0YU4GFpLnM7bJcqABnduSIW+KxsjSXFaWqfLZmD3/5eC37sk9wfa+mPDisNdXjY8MdmjGmHChJc5mvB/+FJ3sNLwMpVfV/5jYTkeruSP+iw86riwiqapNsRQARYfiZDejfMpl/frGBKd9tZeaaPfz5N+0YfmZ96xhgjAkav0kG+MTjfTxwCZBZzLFv4fQkW4qTmDz/Win2XCaiJMXH8tiF7bi0SyMefX81d721jEGtU/jbRe1pUjsh3OEZYyqAEjeXiUgU8JWqDglNSGWnsjaXeVNw8hRTvtvGs1+s56Qq957Vklv7NyM2uqz6hhhjyotQ9y5riTPhpa8A+rozMCMi14nIsyLis4wJr5joKG7pl85XDwxkYKsUnpq5nqHPfsO0xdutF5oxptT8JpnCWZcLX8DHwMN+ir0EHBORjsBDwDbgzdOO1oRcgxpVeeX6bky6sRtJ8bE8+v5qrp24iH1Hjoc7NGNMOeQ3yRTOuuzxaqWq7/kpVuAufXwR8JyqPgck+buWiEwWkX0issZjW20R+VJENro/vS5uLyLnish6EckQkUf8Xcv4dlabenx0d1/+eUVHVu48zNlurcZmdzbGlEQgNZlZgWwrIltEHgWuAz4VkWggkP6xrwPnFtn2CDBLVVsCs9zPReOJBl4AzgPaAleLSNsArmd8EBEu69qYGff2p23D6jz6/mqumvAd6/fYmjXGmMAUm2REJF5EagPJIlLLrVHUFpE0oKGf814FnABuUdU9QCPgaX/BqOpcoGg354uAKe77KThr1BTVA8hQ1c2qmge87ZYzQdAspRrTbu3FU5d1YMPeHIaNnctNry1m28Gj/gsbYyo1X12YbwPG4CSUZR7bj+DUGorlJpZnPT5vB94oZYz1VHW3e57dIlLXyzGNgB0en3cCPUt5PeOFiHBl9yac3bYeUxduY8LczQwbO5ffnt2KG3qnUTUuOtwhGmMiULE1GVV9TlXTgQdVNd3j1VFVx3srU7STQAmXbD4d3kYPen14ICKjRWSJiCzZv39/CEOqmGonxnHPWS354v4B9G2ezBOf/UjfJ2fz0crihk4ZYyqzYmsyIjJEVWcDu0Tk0qL7VfV9L9v8Ptwvhb0i0sCtxTQA9nk5ZifQxONzY4oZMKqqE4AJ4IyTCXawlUWDGlWZNLI73289xOOfruPeacv5fO0eHj3vDBrXsoGcxhiHrwf/A92fv/HyusDfiUWkn4jc5L5PFpH0Usb4EXCj+/5G4EMvx3wPtBSRdBGJA0a45UyIdU+rzbu39+Z3w1rz5dq9DHr6a/4wfTWHjuaFOzRjTAQI1QSZfwa6Aa1VtZWINATeUdW+fspNAwYBycBe4M/AB8B/cQaAbgeuUNVD7jknqupwt+xwYCwQDUxW1cf9xWkj/oMr83AuL329ibcWb6dalRgeOKcV1/RIJcZmDTCmQinJiH+/SUZE7gNew5n08lWgC/CIqn7ho8wKoDOwTFU7u9tWqWqHgH6DMmJJJjQ27M3msY/W8u2mg7RpUJ0/DG9Dv5bJ4Q7LGBMkwZ5W5mZ3ZuVzgLrATcA//JTJcwdjqhtQYiDBmIqhVb0kpo7qyYvXduFIbj7XTVrEdRMXsWGvja8xprIJJMkU9twaDrymqivx3pvL039F5BWgpojcCnyFUwsylUThcgKzHxzIny5oy+pdWZz33Dz++vEPZOXmhzs8Y0wZCaS57DWccSjpQEecZx5fq2pXP+WG4tR+BPhcVb8MSsRBZM1lZefQ0Tye+WI90xZvp0bVWG7t34wb+6RRrUogq00YYyJJsJ/JRAGdgM2qelhE6gCNVHWVl2Nb4AyeXFBk+wBgl6puCvB3KBOWZMreml1ZPPvlBmb/uI8aVWMZ2SeNkX3SqJUYF+7QjDEBCuozGVU9parL3ATzmKoe9JZgXGP531UxAY65+0wl175RDSaP7M6Hd/WlZ3ptnpu1kQFPz+HVuZs5UXAy3OEZY4KspH1LL/SzP81bAlLVJUBaCa9lKrCOTWoy4YZuzBzTn65Na/H4jHWc86+5fLZ6N6HoVm+MCY+SJhl/D/zjfeyrWsJrmUrgjPrVef2mHky5uQdVYqK4Y+oyhj8/nxmrd9uyAsZUACVNMj4f9gPfu73JfkVEbgGWlvBaphIZ2CqFGff259krO3Ii/yR3Tl3GsLFz+XDFLk5asjGm3ArkwX8KcCtOc9fPXYFU9WYvx9YDpgN5/JJUugFxwCXu7MwRwx78R6aTp5RPV+9m3KyNbNyXQ3pyIncOas7FnRsRa7MHGBN2we5d9i0wDydp/Pxk1tfqmCIyGGjvflzrTrTpub+Wqv4USIChZEkmsp06pcxcu4fxszP4YfcRGteqyh2DmnN518ZUibGlBYwJl2AnmRWq2ikYgXmcc5mqdgnmOUvDkkz5oKrM/nEf42ZnsGLHYepXj+fOwc25slsT4mMt2RhT1oI9rcwn7uSTweSvA4ExPxMRzmpTj+l39uHft/SkSe2q/OnDtQx6+msmzttMzomCcIdojClGIDWZbCARZznlfJwEoapavdQXtZqMOQ2qynebDvLcrI0s2nKIpPgYruvVlJv6pFG3uq8OjsaYYChJTcbvnB4hWojMmFITEfq0SKZPi2RW7DjMhLmbeOWbTUyat4VLOjfi5n7ptK5vt60xkcDXyphnqOqPIuK1xqGqy07jutZcZoKiU5OavHhtV7YeOMrE+Zt5Z8lO/rNkB/1bJnNz33QGtkohKspuN2PCpdjmMhGZoKqjRWSOl92qqkN8nlikH9BSVV9zu0FXU9Ut7r7aqnrodIM/XdZcVvEcOprHtMXbeeO7rew9coJmKYnc1Dedy7o0IiHOJuM0JhiC2ruslAGUdmXMyThLO+9T1fbutqdxlnzOAzYBN6nqYS9lt+LMm3YSKAjkC7AkU3HlFZziszW7mTR/C6t2ZlE9Poare6Yysk8aDWrY5BPGnI5ISDIrKMXKmO5szTnAGx5J5hxgtqoWiMiTAKr6sJeyW4Fuqnog0DgtyVR8qsqy7T8xaf4WZq7ZQ3SUcFmXxtw+sDlpybaWnjGlEdQH/6WUp6oqIiVaGVNV54pIWpFtnss8LwQuD1qUpsITEbo2rU3XprXZcegYr87bzNvf7+C/S3ZwfoeG3DmoOW0alLqjpDHGj1DN0RGqlTFvBj4rZp8CX4jIUhEZHYRrmQqmSe0E/npRe+Y/PJhbBzRj9rq9nPfcPEZN+Z5l28M+AYUxFVIg42QEuBZopqp/FZFUoL6qLvZTrlQrY7o1mU8Km8s8tv8B5znPpeolaBFpqKqZIlIX+BK4R1XnejluNDAaIDU1teu2bdsCCctUQFnH8pny3VYmL9jC4WP59G5Wh7sGt6Bvizo4t70xxptgTyvzEnAKGKKqbUSkFvCFqnb3U64pTu+yr0QkAYhWVW8LmhUtl0aRJCMiNwK3A2ep6rEAzvEYkKOqz/g6zp7JGICjJwqYtng7E+ZuZl/2CTo2qcldg5pzdpt61v3ZGC+CPa1MT1W9CzgO4E5s6XOtXLeJ7F3gFXdTI+CDQALycq5zgYeBC4tLMCKSKCJJhe9xalBrSnM9U/kkVolhVP9mzHt4MP/vkjP56Wgeo99cyrnPzeXdpTvJKzgV7hCNKbcCSTL5IhKN88yjcOp/f//V3QX0BY4AqOpGoK6/C4nINOA7oLWI7HTXoRkPJAFfisgKEXnZPbahiMxwi9YD5ovISmAx8KmqzgzgdzPmZ1ViormmZyqzHxjIcyM6IQgPvrOSfk/O5oU5GRw+lhfuEI0pdwJpLrsWuAroAkzB6d31R1V9x0eZRaraU0SWq2pnEYnB6c7sswtzWbPmMuOLqjJ34wEmztvMvI0HqBobzZXdGnNzv3Sa1rHuz6byCloXZhGJArYADwFn4TzEv1hV1/k57zci8nugqtsB4E7g40ACMiZSiAgDW6UwsFUK63YfYeK8Lby1eDtvLNzGsLb1uXVAOl2b1g53mMZEtEBqMt+pau8SndTpmjMKj95lwERvvcLCyWoypqT2HjnOlG+3MnXRdrJy8+mcWpNR/ZoxrF09YmzVTlNJBLt32V+AVcD7gSQJt/azqmgX5EhkScaU1tETBby7dCeTF2xh28FjNKpZlZv6pnFV9yYkxceGOzxjQirYSaZwPZkCnB5mfteTEZGpwKOquj3gqMPAkow5XSdPKV+t28ukeVtYvPUQ1arEMKJ7E0b2TaNxrYRwh2dMSETC3GWzge44Pb2OFm5X1QuDfrHTYEnGBNPKHYeZNH8Ln67ejapy3pkNGNUvnc6ptcIdmjFBFeyazABv272NpvcoM7CYMt8EElRZsSRjQiHzcC5Tvt3KW4u3k328gK5NazGqXzrD2tW3wZ2mQgh2kvHsFRYP9ACW+ltPpjywJGNCKedEAe8s2cHkBVvYcSiX1vWS+O3QlgxrV9+mrTHlWkiby0SkCfCUql7t45hs3MGbHrKAJcADqrq5RBcNEUsypiycPKV8siqT577ayOYDR2nXsDr3D23FkDPqWrIx5VKok4zg9B4708cxfwEygbdwOgqMAOoD64E7VHVQiS4aIpZkTFkqOHmKD1Zk8vysjWw/dIxOTWpy/9BW9G+ZbMnGlCvBbi4bxy+1kiigE7BVVa/zUWaRqvYssm2hqvYSkZWq2jGQ4ELNkowJh/yTp3h36U7GzdpIZtZxuqfV4v6hrendvE64QzMmIMFetMzzr3ABME1VF/gpc0pErsSZJBN+vdBYRA3INKasxUZHcXWPVC7t0oj/fr+D8XMyuPrVhfRpXocHzmllswiYCiWQmsx9qvqcv21F9jcDngN64ySVhcBvgV1AV1Wdf7qBB4PVZEwkOJ5/kqmLtvPS1xkcyMljYKsU7h/aio5NaoY7NGO8CnZz2TJV7VJk23JV7XwaMUYESzImkhzLK+CN77bxyjeb+OlYPn2a1+GWfukMbl3Xuj6biBKUJCMiVwPXAP2AeR67koCTqnq2jwBaAS8B9VS1vYh0wFkP5u8B/g5lwpKMiUQ5Jwr498JtTPl2K7uzjpOenMhNfdO4rEtjEqsE0sJtTGgFK8k0BdKBJ4BHPHZl4/QuK/ARwDfA74BXCms8IrIm0uYzsyRjIln+yVPMXLOHSfO3sGLHYZLiY7i6Ryo39kmjUc2q4Q7PVGKRMK3M96ra3bNZTURWqGqnoF/sNFiSMeXFsu0/MWn+Fmau2QPAue3qc3O/dLqk1rTuz6bMBbV3mYj0AsYBbXCWXY4GjvqaIBM4ICLN+WU1zcuB3YEEZIz5X11Sa9HlmlrsOpzLG99uZdri7Xy6ejcdm9Tkln7pnNe+PrG21ICJQIHcleOBq4GNQFWcdWLG+SlzF/AKcIaI7ALGAHf4u5CITBaRfSKyxmPbYyKyy116eYWIDC+m7Lkisl5EMkTkEW/HGFPeNapZlUeHt+G7R8/ibxe140huPvdOW86Ap+bw0tebbIloE3EC6V22RFW7iciqwuWTReRbVe3j9+QiiUCUqmYHFIwzGWcO8Ebh8xsReQzIUdVnfJSLBjYAQ4GdwPfA1ar6g6/rWXOZKe9OnVK+3rCPSfO3sCDjIIlx0Vzbqymj+qVTt3p8uMMzFVSwB2MeE5E4YIWIPIXT7OV1gXMRub+Y7QCo6rO+LqSqc0UkLYCYiuoBZBTOiSYibwMXAT6TjDHlXVSUMOSMegw5ox7rdh/hlW82MXHeZl7/ditXdG3MbQOak1rH1rUx4RNIc9n17nF346wN0wS4rJhjk9xXN5zmsUbu63ag7WnEebeIrHKb07wtztEI2OHxeae7zZhKo02D6owd0Zk5Dw7isi6NeWfJTgb/82t++58VbNgbUGOCMUHns7nMbYaa4muesmLKfQFcVthMJiJJwDuqem4AZdOATzyay+oBB3A6EfwNaKCqNxcpcwUwTFVHuZ+vB3qo6j1ezj8aGA2Qmpraddu2bSX51YwpN/YeOc7EeZuZumg7x/JOck7betw5uAWdbCYBc5pK0lzmsyajqieBFLe5rCRSAc8nkHlAWgnPURjDXlU9qaqngFdxmsaK2olTwyrUGGcWaG/nm6Cq3VS1W0pKSmlCMqZcqFc9nj+c35YFDw/hvrNasmjLIS5+YQHXTVzEt5sOEIrhC8YUFcgzma3AAhH5iF8vpezr+cqbwGIRmY5TA7kEmFKaAEWkgaoWdn++BFjj5bDvgZYiko4zP9oInNkKjKn0aiXG8duhrbh1QDPeWrSNV+dt4ZpXF9E5tSZ3DmrBWWfYtDUmdALpXfZnb9tV9S9+ynUB+rsf56rqcr/BiEwDBgHJwF7gz+7nTjjJaitwm6ruFpGGwERVHe6WHQ6MxRnHM1lVH/d3PetdZiqj4/kneXfpTl7+ZhM7f3JW7LxzcHPOP7MBMTbWxgQgJCP+RSRRVY/6P7L8sCRjKrOCk6f4eFUmL87ZxMZ9OaTWTuD2gc25rGsjqsREhzs8E8GC9kzGPVlvEfkBWOd+7igiLxZz7LIAzuf3GGNM6MVER3FJ58Z8PmYAr1zflVoJsfx++moGPDWHifM2c/REsdMTGhOwQJrLFuEsOvaRv8kuRSQXZ2aAYk8H1FDV1NKHHDxWkzHmF6rKgoyDvPh1Bt9uOkjNhFhu7pvOjX3SqFE1NtzhmQgS7MGYqOqOIpPwnSzm0DMCOF1xZY0xYSQi9GuZTL+WySzb/hMvzM7g2S838OrczYzsm8bNfdOplVjSjqamsgskyewQkT6Aul2Z78VtOitKVW3QiTEVQJfUWkwa2Z01u7IYPzuDcbMzmDx/C9f3TmNU/3SSq1UJd4imnAikuSwZZynls3Ge4XwO3KeqB0MfXmhZc5kxgVm/J5vxczL4ZFUmVWKiuLZnU24b0MzmR6ukwr6eTHlhScaYktm0P4cX5mTw4YpMoqOEEd2bcPvA5jS0RdQqlWD3LmsmIh+LyH53Gv4PRaRZAOWaisjZ7vuq7tQyxphyrHlKNZ69shOzHxjIpZ0b8dai7Qx8eg6Pvr+KHYeOhTs8E4ECGXn1FvBfoAHQEHgHmOargIjcCryLs6YMONO8fFDqKI0xEaVpnUT+cVkHvnloMCO6p/Le0l0MeuZrHnxnJVsOVKjhdOY0BdSFWVV7Ftm2UFV7+SizAmeOsUUe3Z5Xq+qZpx9y8FhzmTHBsSfrOK/M3cRbi7aTf/IUv+nYkLsHt6BlPWvAqIiC2lwGzBGRR0QkzW0Cewj4VERqi0jtYsqcUNWfJ8gUkRjcpZiNMRVP/Rrx/Pk37Zj/8BBu7d+ML3/Yyzlj53Ln1KX8kHkk3OGZMAqkJrPFx25V1f95PuMubnYYuAG4B7gT+EFV/1D6UIPPajLGhMaho3lMnr+FKd9uJftEAWe3qce9Z7WgQ+Oa4Q7NBEHYe5eJSBRwC3AOzij/z3Ems4yo2owlGWNCKys3n9cXbGXygi1k5eYzqHUK9wxpSdem3tYeNOVFUJOMu3DZ+Tjrwfw8eNPfUsrlgSUZY8pG9vF83ly4jYnztnDoaB59W9ThniEt6dWsTrhDM6UQ7CQzAzgOrAZOFW73NtW/iKzGx7MXVe0QSFBlxZKMMWXrWF4Bby3azitzN7M/+wQ90mpzz1kt6NcimSJTV5kIFuwksyrQ5CAiTX3tj7RpZyzJGBMex/NP8vbi7bz8zWb2HDlOpyY1ufesFgxuXdeSTTkQ7CTzJDBLVb8IRnCRxJKMMeF1osBZQO3FOZvYdTiX9o2qc8+QlgxtU89W64xgwe7CvBCYLiK5InJERLJFxGefRBHpJSLfi0iOiOSJyEl/ZYwxlU+VmGiu7dmUr383iKcu70DO8QJue3Mpw5+fx6erdnPqVET1FTKlEEhNZjNwMbA60N5hIrIEGIEzO0A3nK7MLfx1YRaRycAFwL7C9WpE5D9Aa/eQmsBhVe3kpexWIBtnKYGCQLKs1WSMiSwFJ0/xyardjJ+TQca+HDo0rsHvh7exDgIRJtg1mY3AmpJ2P1bVDCBaVU+q6mvA4ACKvQ6cW+Q8V6lqJzexvAe876P8YPfYgH55Y0xkiYmO4uLOjfh8zACevbIjB7JPMGLCQkZNWULGvpxwh2dKIZD1ZHYDX4vIZ8CJwo1+ujAfc9eeWeEOzNwNJPq7kKrOFZE0b/vEeRp4JTAkgJiNMeVYdJRwaZfGDD+zAZMXbOGlOZsYNnYuI7o3YczZrUhJsvVsyotAajJbgFlAHJDk8fLlevfcdwNHgSbAZaUPE4D+wF5VLW55ZwW+EJGlIjK6uJOIyGgRWSIiS/bv33+aIRljQik+Npo7B7Xg698N4rqeqfzn+x0MenoO42dvJDfPFtktDwIe8S8iiaoa0PSqIpII5KrqKfdzNFBFVf3OBe7WZD4pfCbjsf0lIENV/1lMuYaqmikidYEvgXtUda6va9kzGWPKl837c3hq5npmrt1D/erx3H9OKy7r0pho64lWpoK9nkxvEfkBd8llEekoIi/6KTYLSPD4XBX4KpCAiokhBrgU+E9xx6hqpvtzHzAdZxZoY0wF0iylGi9f35V3bu9N/RrxPPTuKs5/fh5zN1irRKQKpLlsLDAMOAigqiuBAX7KxKvqz0/p3PcJPo7352zgR1Xd6W2niCQWLorm1qLOAdacxvWMMRGse1ptpt/Zh/HXdOZY3klumLyYGyYvZt1uGykRaQJJMqjqjiKb/DWGHhWRLoUfRKQrkOvvOiIyDfgOaC0iO0XkFnfXCIoslCYiDd0pbwDqAfNFZCWwGPhUVWf6u54xpvwSES7o0JAv7x/AH89vw8odhxn+/Dx+985K9mQdD3d4xhXIOJl3gWeB8UAv4F6gm6qO8FGmO/A2kOluagBcpapLgxF0sNgzGWMqjqxj+Yyfs5Ep324DgWt6pHLHoObUqx4f7tAqnGBPK5MMPIfTZCXAF8C9qnrIT7lYnEGUgtPUlR9IQGXJkowxFc+OQ8cYPzuDd5ftJCZKuKZnKncMbE5dSzZBE+wk01dVF/jbVmT/FcBMVc0WkT8CXYC/q+qyQIIqK5ZkjKm4th88xvg5G3lv2S5iooRrezbl9oHNLNkEQbCTzDJV7eJvW5H9q1S1g4j0A54AngF+r6o9AwmqrFiSMabi237wGONmb+T95U6yua5XU24b2Iy6SZZsSqskSabYEf8i0hvoA6SIyP0eu6oD0X7OW9gx4HzgJVX9UEQeCyQgY4wJptQ6CTx9RUfuHtKCcbMzeP3brUxdtI3rejbltoHNbfaAEPPVuywOqIaTiDxH+h8BLvdz3l0i8grONDAzRKSKn2sZY0xINa2TyDNXdGTW/QM5/8yGTF6whf5PzebxT39gf/YJ/ycwpRJIc1nTki42JiIJOBNdrlbVjSLSADgz0takseYyYyqvLQeOMn52BtOX7yQuJorrezk1m+RqVrPxJ6jPZCoySzLGmC0HjjJu9kY+WL6LKjHRXN+7KaMHNLNk44MlmQBZkjHGFNq8P4fxszP4YIWTbG5wk00dSzb/w5JMgCzJGGOK2uQmmw8Lk02fpozub8nGU7C7MMcDtwDtgJ/7/KnqzacTZCSwJGOMKU7GvhzGz97IhyszqRobzQ290xg9oBm1E+PCHVrYBXtlzDeB+jiTZH4DNMZZ5tgYYyqsFnWrMXZEZ7787UCGtq3HK3M30e/J2Tw580cOHc0Ld3jlRiA1meWq2tljgGUs8LmqlvsVKq0mY4wJVMa+bJ6flcHHqzJJiI3mxj5p3Nq/GbUqYc0m2DWZwjnHDotIe6AGkFbK2IwxplxqUTeJ56/uzBdjBjCkTT1e+sap2Tz9+Y/8ZDWbYgVSkxkFvAd0AF7DGaD5J1V9OfThhZbVZIwxpbVhbzbPz9rIp6t3kxAbzci+Ts2mZkLFr9lY77IAWZIxxpyuDXuzeW7WRmas3k1iXAwj+6Qxqn96hU42wV5+uZ6ITBKRz9zPbT0WEzPGmEqtVb0kXrimCzPvG8DAVimMn5NBvyfn8M8v1pN1LOJWOClzgTyTeR34HGjoft4AjAl2ICLSRETmiMg6EVkrIve522uLyJcistH9WauY8ueKyHoRyRCRR4IdnzHG+NK6fhIvXNuFmWP6M6BVMuNmZ3DWs98w58d94Q4trAJJMsmq+l/gFICqFuB/+eXSKAAeUNU2OCtw3iUibYFHgFmq2hKY5X7+FRGJBl4AzgPaAle7ZY0xpkydUb86L17blU/u6UdytThuev17/u+DNeTmheLPZuQLJMkcFZE6gAKISC8gK9iBqOruwkXNVDUbWAc0Ai4CpriHTQEu9lK8B5ChqptVNQ9n6eeLgh2jMcYEqn2jGnxwV19u7Z/Omwu3cf64eazaeTjcYZW5QJLM/cBHQHMRWQC8AdwTyqBEJA3oDCwC6qnqbnASEVDXS5FGwA6PzzvdbcYYEzbxsdH84fy2vDWqJ8dOnOTSF79l/OyNnDxVeTpc+U0ybu1iIM4CZrcB7VR1VagCEpFqOF2mx6jqkUCLednm9V9RREaLyBIRWbJ///7ShmmMMQHr0yKZz8cM4Nz29Xnmiw1c9cp37Dh0LNxhlYlik4yIdBeR+vDzc5iuwOPAP0WkdiiCcWcTeA+Yqqrvu5v3uuvR4P709hRtJ9DE43NjINPbNVR1gqp2U9VuKSkpwQveGGN8qJEQy7irOzP2qk6s35PNuWPn8s6SHVT0YSS+ajKvAHkAIjIA+AdOU1kWMCHYgYiIAJOAdar6rMeuj4Ab3fc3Ah96Kf490FJE0kUkDhjhljPGmIghIlzcuRGfjelP+0Y1+N27q7hz6rIKPWOAryQTraqH3PdXARNU9T1V/T+gRQhi6QtcDwwRkRXuazhOchsqIhuBoe5nRKShiMyAn2tad+N0tV4H/FdV14YgRmOMOW2NayXw1q29eOS8M/hq3V6GjZ3L3A0Vs/m+2BH/IrIG6KSqBSLyIzBaVecW7lPV9mUYZ0jYiH9jTLitzcxizNsr2Lgvh5F90njkvDOIj40Od1g+BWvE/zTgGxH5EMgF5rknb0EIujAbY0xl1K5hDT6+px839U3j9W+38ptx81mzq+L8iS02yajq48ADOCP+++kvVZ4oQtyF2RhjKpP42Gj+/Jt2vHFzD7Jy87nkxQW89PWmCtHV2SbItOYyY0wE+eloHr+fvprP1uyhR3ptnr2yI41rJYQ7rF8J9noyxhhjykitxDhevLYLz1zRkR8yj3De2HlMX76z3HZ1tiRjjDERRkS4vGtjPruvP63rJ/Hb/6zknmnLy+WszpZkjDEmQjWpncB/buvN74a1ZuaaPQwbO5cFGQfCHVaJWJIxxpgIFh0l3DW4BdPv7EtClWiunbiIv33yA8fzy8eszpZkjDGmHDizcQ0+vac/1/dqyqT5W7ho/ALW7Q50esfwsSRjjDHlRNW4aP52cXteG9mdg0fzuGj8AibO28ypCO7qbEnGGGPKmcFn1OXzMf0Z2DqFv3+6jusmLSLzcG64w/LKkowxxpRDdapVYcL1XXnysjNZseMw546dy0crvU4+H1aWZIwxppwSEa7qnsqMe/vTvG417p22nDFvLycrN3K6OluSMcaYci4tOZF3buvNb89uxcerdnPe2Ll8t+lguMMCLMkYY0yFEBMdxX1nt+Td23sTFxPFNRMX8sRn6zhREN6uzpZkjDGmAumcWotP7+3PiO6pvPLNZi5+4Vs27M0OWzyWZIwxpoJJrBLDE5eeyas3dGPfkeNcMG4+k+dvCUtXZ0syxhhTQQ1tW4+ZYwbQr0Uyf/3kB258bTF7so6XaQwRk2REpImIzBGRdSKyVkTuc7c/LSI/isgqEZkuIjWLKb9VRFa7yzbb/P3GGAOkJFVh0o3d+PvF7fl+6yGGjZ3LjNW7y+z6EZNkgALgAVVtA/QC7hKRtsCXQHtV7QBsAB71cY7Bqtop0HUOjDGmIlNVdmflMmf9PrJy82ldvzpZufncOXUZ32zYXyYxxJTJVQKgqruB3e77bBFZBzRS1S88DlsIXB6O+IwxJpKdKDhJxr4cfsg8wrrd2azbfYR1e45w2GN5gNTaCQxrV492DWvQqXHNMokrYpKMJxFJAzoDi4rsuhn4TzHFFPhCRBR4RVUnFHPu0cBogNTU1KDEa4wxZUlVWZt5hG83Hfg5oWTsy6HAfbAfHxtF6/rVOa99fdo0qE6bBtU5o34SSfGxZR5rxCUZEakGvAeMUdUjHtv/gNOkNrWYon1VNVNE6gJfisiPqjq36EFu8pkAzvLLQf8FjDEmRDbtz+HjlZl8tDKTzfuPAlC/ejxtGiQx5Iy6PyeU9OREoqMkzNE6IirJiEgsToKZqqrve2y/EbgAOEuLWYNUVTPdn/tEZDrQA/ifJGOMMeXRmLeX88EKZ26ynum1uaVfOue0rU9KUpUwR+ZbxCQZERFgErBOVZ/12H4u8DAwUFWPFVM2EYhyn+UkAucAfy2DsI0xpkyc274BH6zI5JN7+tG+UY1whxMwKaZiUOZEpB8wD1gNnHI3/x54HqgCFE7Es1BVbxeRhsBEVR0uIs2A6e7+GOAtVX3c3zW7deumS5ZYb2djjCkJEVkaaC/eiKnJqOp8wFsj4oxijs8EhrvvNwMdQxedMcaY0oikcTLGGGMqGEsyxhhjQsaSjDHGmJCxJGOMMSZkLMkYY4wJGUsyxhhjQsaSjDHGmJCJmMGY4SAi+4FtpSiaDBwIcjjBYrGVXKTGBRZbaVlspRNobE1VNSWQE1bqJFNaIrIkUtessdhKLlLjAouttCy20glFbNZcZowxJmQsyRhjjAkZSzKl43VBtAhhsZVcpMYFFltpWWylE/TY7JmMMcaYkLGajDHGmJCp1ElGROJFZLGIrBSRtSLylyL7HxQRFZHkYspvFZHVIrJCRJZ4bK8tIl+KyEb3Z62yjE1EWrsxFb6OiMgYd99jIrLLY9/wYMUW6LlF5FwRWS8iGSLyiMf2kH1vgcQmIk1EZI6IrHPL3uexLxK+tzK/3wL83sJyv7n77nHvpbUi8lQx5cv8fgsktnDdb4HE5h4TnPtNVSvtC2f9mmru+1hgEdDL/dwE+BxnHE1yMeW3etsHPAU84r5/BHiyrGPzOE80sAenXzvAY8CDofjeAjm3G88moBkQB6wE2ob6ewswtgZAF/d9ErDBI7awfm/hut9K+nuX8f02GPgKqOLuqxtB91sgsYXrfvMbWzDvt0pdk1FHjvsx1n0VPqT6F/CQx+eSuAiY4r6fAlwcxtjOAjapamkGnZYmNn96ABmqullV84C3cb4vCP335q/sblVd5r7PBtYBjUoaQyhi8yOs31sRZXm/3QH8Q1VPuMft81I8XPeb39jCeL8F8r35UqLvrVInGQARiRaRFcA+4EtVXSQiFwK7VHWln+IKfCEiS0VktMf2eqq6G5wbCagbhtgKjQCmFdl2t4isEpHJpWkiKC62AM/dCNjh8Xknv/yHFbLvLcDYPM+RBnTG+b+/QuH83iAM91sJYitUlvdbK6C/iCwSkW9EpLuXouG63wKJzfMcaZTd/RZobMG5306nOlaRXkBNYA7QAecfuob6qDK6+xq6P+viVMMHuJ8PFznup7KOzd0fhzNFRD2PbfVwmhCigMeByUGKrX0g5wauACZ6fL4eGBfi7y2g2DzKVQOWApdGyvcWpvutpN9bWd9va4DncZqFegBbcHvMRsD95je2MN5vAcUWrPut0tdkCqnqYeBrnKpgOrBSRLYCjYFlIlLfS5lM9+c+YDrOPxjAXhFpAOD+LGl19LRjc50HLFPVvR7n2quqJ1X1FPCqR8ynG9u5AZ57J84zpUKNgUz3fai+t0BjQ0RigfeAqar6vse5wv29heN+Czg2V5nebzj30vvqWAycwpl7y1NY7rcAYwvL/RZobMG63yp1khGRFBGp6b6vCpwNLFfVuqqapqppOP8gXVR1T5GyiSKSVPgeOAfn/xAAPgJudN/fCHxYlrF5uJoiTReFN4frEo+YTze2HwM89/dASxFJF5E4nOaVj9x9ofreAopNRASYBKxT1WeL7Avr9xam+y3Qf9NCZXq/AR8AQ9ztrfilJuUpLPdbILGF634LMLbg3W+nUwUr7y+c5qflwCr3C/yTl2O24jZJAQ2BGe77ZjhVyJXAWuAPHmXqALOAje7P2mUZm/s5ATiI27Tmsf1NYLV73o+ABsGKrbhze4ltOE5Pmk1l9b0FEhvQD6cdehWwwn0Nj4TvLVz3Wwn+TcNxv8UB/3a3LQOGRND95je2MN5vgcQWtPvNRvwbY4wJmUrdXGaMMSa0LMkYY4wJGUsyxhhjQsaSjDHGmJCxJGOMMSZkLMmYck1EToozS+waEXlHRBKCcM5BItLH4/PtInJDKc+V4/8ov+d4XUQud99PFJG2fo7/OV4RGSkiDUt4vcLZd4O61rs4Mw7nBPu8JrLFhDsAY05Trqp2AhCRqcDtwM8D20QkWlVPlvCcg4Ac4FsAVX05KJEGgaqOCuAYz3hH4oyHyPR+dLEGq2rRgY2nRVUHi8jXwTyniXxWkzEVyTyghVsTmSMib+EMaENEPnAn+lvrOdmfOGuNLBNnzY1Z7kSFtwO/dWtI/cVZ2+NB9/gWIvKVe/wyEWkuItXcssvcGsBFvoJ0R1N/6p5jjYhc5W7/k4h8726b4I4IL1r268KagFsreNw9z0IRqeduf0yc9YYuB7oBU93f5XwRme5xrqEi8n7Ra3i55lYReVKctUkWi0gLd/sVbqwrRWSuu22kiIz3KPuJiAzydw1TcVmSMRWCiMTgzJ212t3UA2eUcmHT0s2q2hXnj+69IlJHRFJw5oW6TFU7Aleo6lbgZeBfqtpJVecVudRU4AX3+D7AbuA4cImqdsFZq+Of3hKEh3OBTFXtqKrtgZnu9vGq2t3dVhW4wM+vnQgsdGOZC9zquVNV3wWWANe6tb0ZQBv39wa4CXjNzzUKHVHVHsB4YKy77U/AMPf6FwZ4HlPJWJIx5V1VcaYyXwJsx5kLCmCxqm7xOO5eEVkJLMSZMLElzgJOcwuPU9VDvi7kzuXUSFWnu8cfV9VjOLPZ/j8RWYWzGFQjnFl0i7MaONutHfRX1Sx3+2Bxpl9fjTO3VDs/v3se8In7fimQ5utgdab3eBO4zp3TqjfwmZ9rFJrm8bO3+34B8LqI3IozY7Ax/8OeyZjy7udnMoXcSsRRj8+DcCYH7K2qx9znAvE4yaEk8yoVVzu5FkgBuqpqvjgzZMcXdxJV3SAiXXHm1HpCRL7AWW3wRaCbqu4Qkcd8ncOVr7/MC3WSwP57fg34GKf29Y6qFgRQBn79Pan7e9wuIj2B84EVItIJKODX//Pq73cwFZzVZExlUANnzYtjInIGTg0G4DtgoIikg7N2ubs9G2c53F9R1SPAThG52D2+itubrQawz00wg4GmvoJxe3sdU9V/A88AXfjlj/EBEakGXF7q3/bXfvW7qDN9eybwR+D1EpznKo+f3wGISHNVXaSqf8KZxbcJzqStnUQkSkSacJpT1Jvyz2oypjKYCdzuNmetx2kyQ1X3u50A3heRKJx1MYbi/J/+u+4D/HuKnOt64BUR+SuQj7Mo1lTgYxFZgjOT7o9+4jkTeFpETrnnuENVD4vIqzhNaVtxpqgPhteBl0UkF6cml+vGm6KqP5TgPFVEZBHO/5he7W57WkRa4tTwZuHM2AvOIlir+WWWX1OJ2SzMxlQybu+v5ao6qZj9W3Ga7Q54+3ya1/4aeFBVl5zuuUz5YM1lxlQiIrIUZ52Rf/s4bD8wS0IwGBNnnZL8YJ7XRDaryRhjjAkZq8kYY4wJGUsyxhhjQsaSjDHGmJCxJGOMMSZkLMkYY4wJGUsyxhhjQub/A8hxm5xmPUIxAAAAAElFTkSuQmCC\n", "text/plain": ["
"]}, "metadata": {"needs_background": "light"}, "output_type": "display_data"}], "source": ["# swap dimensions to temperature data\n", "ocean_data_temp = argo.swap_dims({'n_levels': 'temp'}) \n", "ocean_data_temp.psal.plot(y='temp', yincrease=False, label='Practical salinity')"]}], "metadata": {"celltoolbar": "Tags", "kernelspec": {"display_name": "Python 3", "language": "python", "name": "python3"}, "language_info": {"codemirror_mode": {"name": "ipython", "version": 3}, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.8"}}, "nbformat": 4, "nbformat_minor": 5}